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X-ray radiography is industrially useful in the physical inspection of welded, cast, or as- 
sembled products. X-ray diffraction analysis contributes to the understanding of the details of 
the processes of manufacturing and treatment of alloys, ceramics, glasses, textiles, and other 
materials. Both applications have led to the control and improvement of industrial processes 


and products. 





INTRODUCTION 


HE applications of the physical properties of 

x-rays have led to the development of two 
tools of industrial importance: a tool for seeing 
through opaque objects, and a tool for seeing 
very small things. Of these tools, the first— 
radiography—is probably widely known through 
its medical applications; it is sufficient, therefore, 
to mention only a few examples of its industrial 
uses. On the other hand, the operation of the 
second tool—x-ray diffraction analysis—requires 
some explanation. 


X-RAyYs AS A TOOL FOR SEEING THROUGH 
OPAQUE OBJECTS 


Welded structures 


Probably the most extensive industrial applica- 
tion of the tool for seeing through things has been 
to welded structures. Its advantages are readily 
apparent. It provides an impersonal test with a 
permanent record. This test can be made non- 
destructively on each completed object. Radio- 
graphic inspection, of course, does not attain its 
full significance unless it is related to other 
physical tests. Once this relationship is estab- 


1 Presented at the Conference on Industrial Physics held 
at the University of Pittsburgh on November 15, 1935. 

?The Cooperative X-Ray Laboratory is sponsored 
jointly by the Departments of Physics and of Chemistry 
and Mellon Institute of Industrial Research at the Uni- 
versity of Pittsburgh. 


lished, however, it gives the manufacturer of 
welded structures a basis for confidence in his 
finished products, and allows him to offer guaran- 
tees which cannot be equalled by manufacturers 
using other types of construction. For example, 
no such comprehensive physical inspection on 
each completed object is possible with riveted 
construction. 

These advantages have led the American 
Society of Mechanical Engineers to incorporate 
in their “Rules for the Construction of Unfired 
Pressure Vessels’’ the requirement that each 
Class 1 vessel be inspected radiographically, and 
to make this test optional for other classes of 
vessels. Radiographic inspection has been used 
extensively on the welded penstocks at the 
Boulder Dam, and its use there has been given 
considerable publicity. 

This type of inspection is not limited to the 
factory, however. With modern x-ray equipment 
there should be no serious difficulty in examining 
in this way, for example, every joint in the 
welded framework of a skyscraper. 


Castings 


A second important application is to the 
foundry. Although radiographic inspection has 
been used extensively on aluminum and steel 
castings and forgings, it is not limited to these 
metals. For example, at the University of Pitts- 
burgh successful examinations have been made of 
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red metal or bronze castings containing up to six 
percent lead; of course such inspection requires 
rather penetrating radiation, for lead is com- 
monly used as a partial protection against x-rays. 
The chief usefulness of x-rays in the foundry 
lies in aiding the development of casting pro- 
cedure. They provide a rapid and positive 
method of determining the locations of faults; no 
method of sectioning can be relied upon with 
equal security. On the other hand, this type of 
inspection is usually too costly to be applied to 
every casting, unless its value is high or unless it 
requires expensive machining. Nor should this be 
necessary, for often a relatively few radiographs 
will indicate the way in which to change the posi- 
tions of the gates and risers, or to alter the mold 
in general, so that the number of unsound cast- 
ings can be greatly reduced. Occasional radio- 
graphic inspection should increase the manu- 
facturer’s confidence in his products and provide 
a way to convince purchasers of their soundness. 


Assembled articles 


Another application is to the inspection of 
assembled articles in which the parts are hidden, 
such as metallic radio tubes and embedded 
electric heating elements. Airplane propellers 
have been inspected in this way. Other applica- 
tions will readily suggest themselves. 


Visual examination 


It is not always necessary to prepare a radio- 
graph. The x-rays transmitted by the object 
under inspection may be made to illuminate a 
sensitive screen, so that the faults may be sought 
visually. This offers a rapid and cheap method of 
examination of many types of completed prod- 
ucts. Such inspection has been adopted by cer- 
tain manufacturers of chewing tobacco, arc-light 
carbons, and some food products. Protection is 
thereby offered, not only to the customer, but 
also to the manufacturer. 

It may be argued that such inspection is not 
warranted by the number of articles rejected. 
Although this objection may be well taken from 
an engineering standpoint, the advertising value 
of such inspection may often be sufficient to 
justify its cost. For example, advertisements of 
national scope offer a radiograph with each 
tobacco pipe of a well-known brand. 
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X-Rays AS A TOOL FOR SEEING SMALL THINGS 


Although the tool, provided by x-rays, for 
seeing very small things is probably not so well 
known as the tool for seeing through things, its 
applications are certainly equally important. 


Physical principles 


The physical basis on which this tool rests 
can be described briefly. Many physical studies 
have shown that the atoms of most substances 
are spherical or approximately so. One is not far 
wrong, then, in thinking of a solid object as 
composed of atoms packed together as though 
they were small rubber balls of different sizes. 
The way in which the atoms of some typical 
substance pack together in the solid form is 
illustrated by Fig. 1. For the present purposes, 
the details of the arrangements are not impor- 
tant; however, it should be observed that a 
regularity of structure exists, and that this 
regularity may extend over thousands to millions 
or more of atoms in one solid body. As a conse- 
quence of this regularity, the atoms lie in 
uniformly spaced parallel planes. In the model of 
rocksalt, for example, planes which contain only 
sodium atoms may be readily selected, and 
parallel to these planes are found others which 
contain only chlorine atoms; others contain both 
types of atoms, such as the planes through the 
diagonals of the square faces. These planes are in 
the neighborhood of one-two hundred millionth 
of an inch apart. 

It is known that x-rays can be reflected from 
these planes just as light is reflected from a 
silvered surface, but there is one aspect in which 
the reflection differs. Because of the penetration 
of x-rays into solids, reflections take place not 
only from one plane, but from large numbers of 
parallel planes. It follows, as a consequence, that 
X-rays cannot be reflected from the atomic 
planes in all directions, as can light from a mirror, 
but reflections occur only when the beam strikes 
these planes at certain sharply defined angles. It 
is only necessary, therefore, to measure the 
angles at which x-rays of known wave-length may 
be reflected in order to determine the distance 
between the atomic planes, and thus to learn 
something, at least, of the atomic arrangements. 

A convenient method by which the reflections 
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Fic. 1. Some of the more common types of crystal lattices. 


(or more strictly, the diffractions) of x-rays at 
angles of incidence from approximately zero to 
ninety degrees on atomic planes are recorded 
simultaneously on a photographic film is known 
as the Debye-Scherrer-Hull method; a descrip- 
tion of it can be found elsewhere.’ The film, when 
finished, is known as a diffraction pattern. or 
photogram. For the present purposes it is 
sufficient to observe that a change in the positions 
of the reflections without a change in the relative 
configuration of the diffraction patterns indicates 
an alteration in only the dimensions of the struc- 
ture, the general arrangement being preserved; 
the appearance of a new configuration of reflec- 


> Cf. George L. Clark, Applied X-Rays, second edition, 
McGraw-Hill Book Company, New York (1932), p. 199. 


tions indicates a new arrangement of atoms, and 
thus, usually, gives evidence of a new substance. 


Metallurgical applications 

This method has had a wide application in the 
field of metallurgy; an example of its application 
to the silver-cadmium alloy system by Astrand 
and Westgren‘ will serve to illustrate the explana- 
tion just given. The phase diagram, shown in 
Fig. 2, is characterized by a large field of solid 
solution of cadmium in silver, a narrow field of a 
second phase which overlaps two others, a third 
and a fourth phase both of which have a range of 
solid solubility for cadmium of some extent, and 
finally a solid solution of silver in cadmium. 


*H. Astrand and A. Westgren, Zeits. f. anorg. allgem. 
Chemie 175, 90 (1928). 
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Fic. 2. The silver-cadmium alloy system. 


The spatial arrangements of the atoms deduced 
from diffraction patterns of these alloys are 
shewn in Fig. 2, below the phase diagram. Over 
the first large range of solid solution, the original 
arrangement of the silver atoms is maintained, 
cadmium atoms being substituted at random for 
the silver atoms, with a resulting expansion of the 
atomic arrangement. In the second phase, the 
cadmium atoms have definite positions with 
respect to the silver atoms, and the binary com- 
pound AgCd is indicated. The third phase has a 
large, complicated structure requiring 52 atoms 
to represent a single structural unit, and foreign 
atoms can find room in it with only a small 
change in the arrangement. In the ideal case this 
structure represents an intermetallic compound 
of the formula Ag;Cds. The fourth phase is 
similar in structure to the last one, cadmium. The 
structure of the fourth phase, in the ideal case, 
also represents a compound which has the formula 
AgCd;. The large structural unit of the third 
phase is similar to that of y-brass, and, in fact, 
this discussion applies just as well to the copper- 
zinc alloy system insofar as the general results are 
concerned. 

Such studies as this allow one to determine 
the type of solid solution, that is, to distinguish 
between the cases in which the solute atoms 


substitute for the solvent atoms and those in 
which the solute atoms slip in between the 
solvent atoms. It also permits the determination 
of the range of solid solution. Further, it gives 
positive proof of the existence of intermetallic 
compounds. Important as these contributions 
are, however, this method still has one greater 
advantage, for it makes possible a classification 
of alloy systems on a new basis. The binary 
alloys of copper, silver, or gold with aluminum, 
zinc, cadmium, or tin have one or more phases, 
the structures of which are similar to those of the 
silver-cadmium system; other systems also fall 
into this classification. When the initial com- 
ponents have similar structures, they are gener- 
ally miscible in the solid state in all proportions, 
but if the structures are different, they are 
miscible, without the formation of compounds or 
simple mixtures, over only a limited range. It is 
the atomic arrangements of the phases, of course, 
which determine the physical properties of the 
alloys. On the basis of information now available, 
it is possible, in many cases, to predict the atomic 
arrangement and the type of phase diagram to be 
expected from a knowledge of the elements 
forming the alloy, and thus to determine in 
advance some indication of the physical prop- 
erties of a new alloy. 

X-ray diffraction methods also give informa- 
tion concerning the processes of rolling and 
drawing. It has been found that there is a close 
relation between the alignment or orientation of 
the crystallites produced in a metal by such 
processes and the resultant mechanical and 
magnetic properties. X-ray diffraction methods 
permit the measurement of the degree of crystal- 
lite orientation. Fig. 3A shows a diffraction 
pattern of an iron sheet with a small amount of 
reduction in thickness produced by cold-rolling; 
the pattern is characterized by nearly continuous 
rings which indicate a random crystallite orienta- 
tion. A diffraction pattern of a similar sheet 
which has undergone a severe reduction is shown 
in Fig. 3B; the rings are broken into short arcs, 
the lengths of which are related to the degree of 
perfection or orientation of the crystallites. This 
material has physical properties which vary with 
the direction in the rolled sheet. It also is possible 
from patterns of this type to estimate the internal 
stresses in certain materials. Such studies as these 
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Fic. 3. Diffraction patterns of ro'led steel sheets. A, low reduction, B, high reduction. 


have led to the improvement of transformer 
steels, and to the development of better rolling 
technique. 

Annealing of rolled or drawn materials at the 
proper temperatures destroys this orientation 
and increases the size of the crystallites; along 
with these changes the properties of the metal are 
altered. It is possible to follow the growth and 
reorientation of the crystallites by x-ray diffrac- 
tion methods. 

Again, x-ray diffraction analysis may be used 
to determine the size of crystallites when these 
are below the range of microscopic visibility, as 
the widths of the reflections are related to the 
sizes of the reflecting crystallites. 


Applications to ceramics 


Analogous applications are found in the field of 
ceramics. Here, however, one deals with systems 
the components of which are compounds instead 
of elements as in the field of alloys. This presents 
no difficulty, however. As an example, it is 
interesting to observe that the components of 
clays may often be identified by x-ray diffraction 
analysis when microscopic methods prove very 
difficult. 


Applications to studies of glasses 


X-ray diffraction analysis has made notable 
contributions to the glass industry. By this 
means it has been shown that the coloring agents 
in some glasses are actually crystalline, even 
though they are very small. Certain opacifiers in 
glasses and enamels are also found to be crystal- 
line. Hence the primary causes of color and 
Opacity may be thought to be associated with 
extremely small crystallites of the proper sub- 


stances. Thus, there is offered a means of under- 
standing and controlling color and opacity. In a 
similar manner, x-ray diffraction analysis is 
pointing the way to the understanding of the 
physical nature of a clear glass. 


Applications to studies of textiles 


Textile manufacturers are also concerned with 
problems of crystallite orientation. The per- 
fection of orientation of the crystallites of fibers 
may be measured in the same way as for rolled 
metals. The diffraction pattern of an unstretched 
cotton fiber, for example, is not widely different, 
in the general character of the reflections, from 
that of lightly rolled steel. Stretching the fibers 
increases the perfection of orientation. The dif- 
fraction pattern of ramie, a long-fibered cotton, 
shows relatively sharp spots, indicating a high 
degree of crystallite orientation. Similar studies 
have shown the way in which mercerization and 
dyeing take place. It has also been shown that 
water acts upon the individual molecules of wool, 
while in cotton fibers it remains outside the 
boundaries of large groups of molecules; this has 
led to an understanding of the cause of wool 
shrinkage, and its reduction in manufacturing 
processes. In a similar manner the steps involved 
in the spinning of artificial fibers may be fol- 
lowed; x-ray studies have led to the betterment 
of spinning processes and the improvement of the 
quality of the spun materials. 


Studies of industrial dusts 


Finally, x-ray diffraction analysis may be used 
as a method of identifying compounds, just as the 
spectroscope is used to identify elements. Fig. 4 
shows an example obtained in the study of the 
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Quartz 


! Factory Dust 


Fic. 4. Diffraction patterns obtained in a study of industrial dusts. 


industrial disease known as silicosis, for which 
quartz is often held responsible. The central 
pattern is that of quartz, below it is the pattern 
obtained from a factory dust, and above it is the 
pattern of the material removed, without heat or 
acid treatment, from a lung of a pneumonia 
patient after autopsy. From the similarity of the 
positions and intensities of the lines, the con- 
clusion is obvious. 


CONCLUDING REMARKS 


X-rays, used as a tool for seeing through 
things, have their industrial applications pri- 


’The sample used in obtaining this photogram was 
prepared by Dr. Olive McElroy at the University of 
Pittsburgh. 


marily in the inspection of fabricated materials. 
X-rays, as a tool for seeing small things, have 
their applications to the understanding of indus- 
trial processes. Both tools have led to the control 
of industrial processes, to the reduction of indus- 
trial wastes, to the solution of industrial research 
problems, to the improvement of products, and 
occasionally to the development of new products. 
Many applications could be suggested besides 
those given. X-rays will not solve every problem, 
nor is such power claimed for them. It is prob- 
able, however, that the applications of x-ray 
methods have succeeded many more times than 
they have failed. Often valuable results have 
been obtained, results which were not even 
suspected at the beginning of an investigation. 
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Some Applications of Physics to Air Conditioning 


EMERSON M. PuGu, Carnegie Institute of Technology 
(Received December 13, 1935) 


HE term air conditioning is applied to the 
cleaning of air and to the control of its 

humidity and temperature. Physics has con- 
tributed much to each of these branches, but I 
will confine my remarks to temperature control 
because of my greater interest in that phase of 
the subject. j 

Some time ago friends of mine purchased an 
electric refrigerator. The increase in their electric 
bill was several times as large as the salesman 
had led them to expect. One day they noticed 
that the maid was leaving the door of the re- 
frigerator wide open most of the time. When 
questioned, she said she was just trying to keep 
the kitchen cool. It is surprising to find how 
many people, even engineers and scientists, 
when they first hear this story do not realize 
that such procedure heats rather than cools the 
kitchen. Assume, for example, a _ refrigerator 
whose electrical power consumption under these 
conditions amounts in heat units to 2 B.T.U. 
per minute, and whose evaporator coil inside 
the box absorbs 6 B.T.U. per minute. The con- 
denser coils must then give off 8 B.T.U. per 
minute to the kitchen. The net result is the 
addition of 2 B.T.U. per minute to the kitchen. 

This simple example illustrates how a mechan- 
ical refrigerator (which is a heat engine running 
backwards) acts as a “heat pump.” Put the 
evaporator in a building and the condenser 
outside and we have a cooling system for summer 
use. Interchange the condenser and the evapo- 
rator and we have a heating unit which will 
pump heat into the building from a colder 
exterior. Our simple illustration shows that the 
heat delivered to the building may be several 
times as great as the heat delivered by the usual 
electrical resistance heater. The efficiency of the 
“heat pump” for heating purposes is the recip- 
rocal of the heat engine efficiency or 7; (7,—T>). 
These facts were very clearly pointed out by 
Kelvin in 1852.! 


Wm. Thomson, ‘‘On the Economy of the Heating or 
Cooling of Buildings by Means of Currents of Air,” 
Glasgow Phil. Soc. Proc. 3, 269 (1852). 
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The advantage to be gained by the great 
efficiency of the ‘‘heat pump”’ for heating build- 
ings has been counterbalanced, until recently, 
by the high cost and size of the equipment 
needed. However, the same equipment can be 
used for both heating and cooling. If we already 
have equipment for cooling, there are many cases 
where it is advisable to use it for heating. This 
is especially true in climates where the winter 
temperatures are not too low, or where electric 
power can be obtained cheaply. 

The “‘heat pump”’ principle has been tried out 
in enough places to prove its usefulness under 
certain conditions. The California Edison Com- 
pany has an office building in Los Angeles which 
has been successfully heated and cooled by this 
principle for some time. The American Gas and 
Electric Company has an office building in 
Salem, New Jersey, which has been successfully 
heated and cooled for more than a year. They 
report? that the heat delivered during the last 
winter season averaged nearly four times the 
electrical input to the system. Instead of using 
the outdoor atmosphere as a heat source in 
winter and a heat sink in summer, this installa- 
tion uses water pumped from a deep well. The 
turning of a few valves changes the summer 
cooling system to a winter heating system. The 
Baldwin-Southwark Corporation’ now makes a 
small compact system suitable for controlling 
temperatures in small rooms. It has operated 
successfully on an electric railway coach running 
out of Philadelphia. This system can be con- 
trolled by a thermostat so that it will either heat 
or cool the coach as required. 

All of the “heat pump” systems which have 
come to my notice so far are using electrical 
energy to operate their compressors. It seems to 
me that some cheaper source of power might 
have great advantages particularly for heating. 


2 Sporn and McLenegan, ‘‘All Electric Heating, Cooling 
and Air Conditioning System,” Heating, Piping and Air 
Conditioning 7, 402 (1935). 

3 Galson, ‘‘Using the Reversed Cycle Refrigerating Prin- 
ciple for a Self-Contained Heating and Cooling Unit,” 
Heating, Piping and Air Conditioning 7, 497 (1935). 
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Fic. 1. Ratio of the amplitude of the internal surface temperature 6; to the amplitude of the 
external surface temperature 9, for sinusoidal temperature variations. 


By burning gas or fuel oil in an internal com- 
bustion engine operating the compressor, more 
heat could be delivered to the building than the 
total heat content of the fuel burned. This is 
true because all of the energy losses in such an 
engine appear as heat which could also be 
utilized for heating. The principle of the gas 
refrigerator might also be used for heating and 
cooling buildings. The fact that such a re- 
frigerator has no moving parts might be a distinct 
advantage. 

Unfortunately the heating efficiency of the 
“heat pump” system decreases as the outside 
use is, therefore, 


temperature decreases. Its 


limited at present to the milder climates or to 
locations where plenty of water is available. 
Because summer heat travels in waves into the 
earth some localities have the ideal situation of 
having well water which is warmest in midwinter 
and coldest in midsummer. Such places may be 
ideal for heat pump installations. 

In the design of cooling systems, the nature of 
heat flow into a building needs detailed c¢on- 
sideration. While the average flow of heat 
through a wall can be calculated in the usual 
manner from the average internal and external 
temperatures, it should be pointed out that this 
is not the significant quantity. What the designer 
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Fic. 2. Time lag to of the maximum internal surface temperature behind the maximum external 
surface temperature, for sinusoidal temperature variations. w=22/7, where T is the period. 


needs to know is the maximum heat flow which 
cannot be calculated so simply. 

Not long ago, a Milwaukee company built a 
movie theater and made the southern wali four 
bricks thick because it was exposed to the sun 
all day. Too late they discovered that this thick 
wall retarded the heat wave so that the peak 
arrived in the building in the early evening when 
the theatre was crowded. Since there was no 
afternoon show, a thin wall of sheet iron would 
have been more effective. 

The heat flow due to the sun effect through 


different types of roof panels has been measured‘ 
by the Research Laboratory of the American 
Society of Heating and Ventilating Engineers in 
Pittsburgh. My mathematical solution of the 
problem, amply verified by these experiments, 
has been published in two journals.‘:® The 
solution seemed too complex for the average 


* Houghton, Blackshaw, Pugh and McDermott, ‘‘Heat 
Transmission as Influenced by Heat Capacity and Solar 
Radiation,” Trans. Am. Soc. Heat. and Vent. Eng. 38, 
231 (1932). 

5 Pugh and Houghton, ‘‘Heat Waves through Insulating 
Walls,” Physics 2, 93 (1932). 
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Fic. 3. Typical curve of the variation with time of the temperature of the surface of roof 
panels exposed to the summer sun. Dotted line shows the sum of two sine curves whose fre- 
quencies are in the ratio of 1 to 2 and whose amplitudes are in the ratio of 1 to 0.32, respectively. 


engineer to use, and so a simpler approximate 
solution was also given in these articles. However, 
the approximate solution was still so complex 
that most engineers continued to use “rule of 
thumb” methods found in the handbooks, even 
though these methods might lead to the kind of 
difficulties experienced by the Milwaukee con- 
cern. Fortunately, since this material was pub- 
lished, it has been found possible to reduce the 
accurate solution to a series of curves so that 
practical problems can now be solved with little 
difficulty. 

We can divide the total heat flow through a 
wall or roof, into two components; one steady, 
which is that usually calculated; the other 
periodic, requiring special treatment. Over a 
twenty-four hour period, the average heat flow 
into a building due to the periodic component is 
zero. The conventional calculation, therefore, 
does give the correct average heat load. However, 
in the design of cooling equipment the peak of 
the heat flow is the important factor and this 
may be several times the average value. It is 
possible to design a wall so that the cooling 
portion of the heat wave coincides with the peak 
of the heat load from other sources, thus reducing 
the cooling requirements. 

If the periodic temperature component of the 
outside wall surface is assumed to be = 45 cos wrt, 
then the periodic temperature component of the 
inside wall surface is given by 


6=0,{ De®! cos (wt+a—~ i) 


+ Ee-*' cos (wi+é—Bl)}, (1) 


where / is the wall thickness. The constants D, 
E, 8, wand 6 depend upon the physical properties 
of the wall (conductivity K, density p, specific 
heat c, inner film conductance e) and the period 
P=27/w. Eq. (1) may be written 


§=80; cos w(t—to), 


where 6; and f) are complex functions of D, E 
B, a, 6,1, w and 6. If a new constant f(=e/K@£) 
is introduced, 6;/8) and t) can be expressed in 
terms of the two dimensionless quantities f and 
Bl, where B=(cpw/2K)'. It will be noticed that 
6; 0) is the ratio of the amplitudes of the internal 
and external temperatures and f) is the time lag 
of the maximum internal temperature behind the 
maximum external temperature. In Fig. 1 4;/4 
has been plotted for various values of f and (i. 
Thus, if the amplitude 4, of the outside tempera- 
ture and the constants f and @/ of our particular 
wall are known, we can find 6; directly from the 
curve. The maximum heat flow is then ¢6; plus 
the average heat flow calculated in the usual 
manner. The time of arrival of this maximum 
heat flow can be obtained from Fig. 2 which 
gives the lag of the inside maximum temperature 
behind that of the outside maximum tempera- 
ture. In order that the curve may be useful for 
all possible frequencies of the impressed heat 
wave, wf) has been plotted instead of the time 
lag to. w= 22/T, where T is the period. 

Figs. 1 and 2 give all the information needed 
for sinusoidal temperatures, but the sun does 
not produce simple sinusoidal temperatures on 
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TABLE I. Physical properties of some panels tested by the Research Bureau of the 
American Society of Heating and Ventilating Engineers, Pittsburgh, Pa. 
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walls or roofs. The actual temperatures follow 
complex curves with periods of twenty-four 
hours. These curves can be analyzed into their 
sinusoidal components which must be treated 
separately, since the values of f and £/ will be 
different for each component. The experiments 
of the Research Bureau under F. C. Houghton‘ 
showed that for practical purposes the tempera- 
tures of surfaces exposed to the sun can be 
expressed by a Fourier series with only two 
periodic terms, the fundamental and the first 
harmonic (double frequency term). The analysis 
of the experimental curvest shows them to be 
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surprisingly similar even when they were taken 
at different times of the year. In particular, the 
percentage of the first harmonic required to fit 
the curves of external temperatures is quite 
constant and averages 32 percent of the funda- 
mental. Also if the peak of the fundamental is 
made to lead the peak of the first harmonic by 
45 minutes (time) the resulting composite curve 
represents the typical experimental curves quite 
well. See Fig. 3. 

If we assume that the composite curve is the 
typical external temperature curve for the sun 
effect, we find we can express everything in 
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Fic. 5. Time lag of the maximum internal temperature 


9m to the range of temperatures of the external surface Ro, 
where external temperatures are assumed to follow curves 
similar to the dotted curve in Fig. 3. 


behind the maximum external temperature, where external 
temperatures are assumed to follow curves similar to the 
dotted curve in Fig. 3. 
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terms of the f and §/ of the fundamental and the 
range of the external temperatures Ro. To calcu- 
late the maximum flow of heat we need the 
maximum of the periodic temperature component 
of the inside surface, @,,. The maximum heat flow 
is given by «6, added to the average heat flow. 
Fig. 4 shows the ratio of 0, to Ro plotted as 
functions of f and @/. Fig. 5 shows the time lag 
of @, behind the external peak temperature. 
Table I gives values of the constants for some 
of the roof panels tested by the Research 
Bureau.‘ Notice that for the 8-inch cork-iron 
panel f=1.65 and 6/=5.1 giving a time lag of 
18.5 hours. 

While the equations given for the constants 
f and #l hold in any consistent system of units, 
they must be modified for the British Engi- 
neering system which mixes feet and inches 
in its physical constants. When K is in B.T.U. 
inch/ft.2-hr.-F°, p is in lbs./ft.*, € is in B.T.U./ 
ft.2-hr.-F°, and / is in inches; the constants be- 
come 8=((2/288)(cp/K))'=0.1044(cp/K)! and 
f=e/KB for a fundamental period of 24 hours. 

It must be remembered that while the curves 
of Figs. 1 and 2 are exact calculations, the curves 
in Figs. 4 and 5 are based upon experiments 
performed in Pittsburgh, Pennsylvania on roof 
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panels. It is not known, therefore, how closely 
the last two curves apply to roofs in other 
latitudes or to side walls anywhere. Neverthe- 
less, it is believed that if they are used on any 
wall or roof in the temperate zones the results 
will be better than can be obtained from present 
handbook methods. Wherever the sun effect is 
important the surface temperatures must follow 
curves somewhat similar to Fig. 3, with a high 
peak in the daytime and a flat cooling curve at 
night. 

When a wall contains many windows the direct 
rays of the sun may heat up the interior so that 
the heat wave through the wall may be largely 
masked. The wave exists, however, and has an 
additive effect. 

It is possible to extend this work to give 
solutions for compound panels made up of slabs 
of different materials. The time required for 
preparing curves for this purpose has not yet 
been available. 

I wish to acknowledge my indebtedness to the 
many who have aided me in this work. I am 
particularly indebted to Messrs. F. O. Mortlock, 
I. Gazdik, and J. J. Dyke. Mr. Mortlock and 
Mr. Gazdik performed many of the calculations 
and Mr. Dyke performed them all independently. 
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How Physics Is Applied in the Oil Industry* 


Pau D. Foote, Gulf Research & Development Corporation, Pittsburgh, Pennsyivania 
(Received December 14, 1935) 


SHALL try to present a bird’s-eye view of the 

general types of problem with which the 
physicist must deal in the various branches of 
the petroleum industry, including discovery of 
oil fields, drilling, production, transportation, 
refining, marketing, and studies of the physical 
properties of the finished products. 

A typical oil reservoir consists of porous and 
permeable sedimentary rock strata originally 
laid down as sand or calcareous deposits in the 
relatively shallow waters which once covered the 
geological basins now forming the petroleum 
producing areas. Through geologic time these 
sediments have been covered by thousands of 
feet of other sediments and the resultant pres- 
sure together with the precipitation of salts from 
the underground waters circulating through 
them has resulted in the consolidation of the 
sediments and their transformation into the 
rocks,constituting the sedimentary land masses. 
Coincident with this gradual transformation has 
been the distortion of these strata into folds and 
dome-like warpings resembling buried hills 
forming the structural features of a typical 
petroleum reservoir. Along with this great bulk 
of material were deposited certain sediments 
containing the remains of marine fauna which 
have been transformed into petroleum hydro- 
carbons by processes not yet understood. 

Since some strata are sufficiently permeable to 
permit the flow of entrained water, this has been 
free to migrate, in compliance with pressure con- 
ditions imposed upon it, and accompanying the 
water has been the oil and gas forced by com- 
paction from the shales into the more permeable 
sands and limestones. Thus one may postulate a 
sheet of water migrating through geologic time 
with an almost infinitesimally low velocity. 
Eventually the fluids reached one of the dome- 
like folds and there, because of the lower density 
of the oil and gas compared with water, the 
former will be trapped. Any excess of gas not 


_* Presented at the Conference on Industrial Physics, 
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soluble in the oil at the existing pressure and 
temperature conditions will segregate as free gas 
in the uppermost portion of the structure and 
the oil will occupy the intermediate zone be- 
tween the gas cap and the water. 

It is the problem of the geophysicist to find 
these buried structures at depths ranging from a 
few thousand feet to two miles or more, when 
there is no suspicion of their presence from sur- 
face indications. Any conceivable physical prop- 
erty of the rocks of the earth's crust that may be 
measured from the surface through thousands of 
feet of cover might be made the basis for a 
method of geophysical prospecting. Of the many 
schemes devised by the physicist, three funda- 
mental methods have met with especial suc- 
cess. These are the gravitational method which 
depends upon differences in rock density; the 
magnetic method which depends upon differ- 
ences in the magnetization of rocks; and the 
seismic method which depends upon differences 
in the rate of propagation of a compressional 
wave or upon the reflection of such a wave at a 
change in lithologic character of the rock. The 
physicist must devise practical methods and 
instruments for accurately measuring these 
effects at the surface and he must develop the 
technique for rapidly and economically in- 
terpreting the field data in terms of geologic 
structure. 

This latter part of the work is no small prob- 
lem. Newtonian potential function theory is of 
course the basis of the interpretation of gravity 
and magnetic data. However, one immediately 
becomes involved with analytical expressions 
which require a prohibitive amount of time for 
computation. Methods for handling the data in 
a direct manner must be developed and machines 
must be devised which can be used to solve 
equations quickly. For example, one integrating 
machine we have constructed does the work of at 
least 350 computers working with ordinary 
mathematical tables. Special transparent com- 
puting charts and other short cut schemes, some 
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involving integral equations, series expansions 
in very peculiar forms and arrays, and methods 
for adapting the ordinary commercial comput- 
ing machines are necessary in order that the work 
may be made effective at reasonable cost. 

The gravity work is conducted with the torsion 
balance, the pendulum, and the gravimeter. 
The torsion balance is the only instrument we 
use which has not been entirely developed by our 
physicists and constructed in our instrument 
shop. This instrument measures the second 
differentials of the Newtonian potential function, 
the two most useful derived quantities being 
curvature and the horizontal gradient of gravity. 
The instrument reads directly to 10-° dynes/g 
cm which is an extremely minute quantity. 
Many problems naturally arise in securing 
measurements to this degree of refinement. One 
incidental problem we solved was the manu- 
facture of suspension wires having a zero tem- 
perature coefficient of zero drift. 

We are using half-second pendulums to de- 
termine gravity to a precision of one part in ten 
million. The pendulums, two to a case swinging 
at 180° phase angle, in vacuum, are constructed 
of quartz with aluminized mirror surfaces. A 
field determination of this precision is made in 
the course of 1/2 hour. Time must be measured 
to one part in twenty million, which is far beyond 
the accuracy of the Naval Observatory signals. 
This is accomplished by using a set of our pendu- 
lums as the standard clock and broadcasting 
time from one of our dozen radio base stations. 
While we do not know the absolute accuracy of 
our time scale, we do have a constant and 
reproducible time scale of the required precision. 
Making gravity measurements to 10~’ is a feat 
for physicists under the favorable conditions 
of the laboratory. It must always be remembered 
that these as well as all geophysical data must be 
obtained with rugged instruments capable of 
being readily transported by auto truck, by mule 
pack, and by manual labor, through swamps and 
across prairie. There are few roads where geo- 
physical prospecting is conducted; these come 
later, after the drilling of the first well. The 
condition of portability of geophysical equipment 
materially controls the design of apparatus and 
taxes the ingenuity of the physicist to the 
highest degree. In passing it is of interest to 
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note that gravity measurements to one part in 
ten million require that the radius of curvature 
of the knife-edge remain constant to about 
100 molecules. The effect of knife-edge wear is 
reduced to a minimum by designing the pendu- 
lum in a peculiar shape such that the radius of 
gyration about the center of gravity is equal to 
the distance between the knife-edge and the 
center of gravity. Also solutions of the differen- 
tial equations of motion show that the mean 
period of a synchronous system of two pendu- 
lums swinging 180° apart in phase, is independent 
of horizontal seismic disturbances. These two 
factors in design have been responsible for a 
precision increase of twenty-fold. 

The gravimeter is simply the familiar fish- 
scale balance, a weight hung on the end of a coil 
spring. When gravity increases at a certain 
location, the weight is deflected a larger amount 
so that the magnitude of the deflectional 
change is a measure of the gravity anomaly. But 
to detect a gravity difference of one part in ten 
million or better with any reasonable size of 
apparatus, requires the measurement of a change 
in deflection of the weight which is equivalent to 
about one-one hundredth the resolving power of 
the best microscope. The development of a 
practical instrument of this type is a real prob- 
lem in instrumental physics. With our gravity 
field instruments we are able to measure the 
daily variation of gravity due to the position of 
the sun and moon and to interpret these data in 
terms of tidal forces of the earth. At Pittsburgh, 
for example, the daily land tides amount to 
from 12 to 23 inches depending upon the phase 
positions of the sun and moon, and there is a 
lag of the actual tide behind the sun and moon 
position, indicating a plasticity of the earth’s 
crust. 

The above-described precision of gravity in- 
struments is necessary when it is recalled that a 
geologic structure giving rise to a major oi! 
field may involve a total gravity anomaly of 
from 0.001 to 0.01 dyne/g or 1 to 10 parts per 
million of gravity. 

In the magnetic method of geophysical pros- 
pecting we measure changes in the vertical com- 
ponent of the earth’s magnetic field with a 
precision of 10-° gauss. This precision requires 
the development of special apparatus of the 
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variometer type. The interpretation of the final 
data is further complicated in that we have the 
vector quantity, polarization, with which to 
deal, and must distinguish changes in magnetic 
intensity due to geologic structure and changes 
due to variations in magnetic susceptibility. 
These latter effects are generally much larger 
than the structure effect which is the objective 
of the survey. 

Two basic methods, with numerous variations, 
are used in seismic prospecting. The most 
generally employed method is reflection shooting, 
the counterpart of the echo method for deep 
sea sounding. With the ocean bottom, a single 
reflection is assured, while over land, reflection 
occurs at every lithologic horizon. Knowing the 
average vertical velocity of sound waves in a 
given area and observing the echo-time-interval, 
the depth to any particular reflecting horizon 
may be determined. This method is particularly 
“useful in areas where lithologic contrasts are 
pronounced, that is, where thick beds of lime- 
stone alternate with shale. The reflection seismo- 
graph record provides only circumstantial evi- 
dence regarding the kinds of rock responsible for 
the reflections. Quantitative information con- 
cerning the properties of buried rocks is obtained 
by refraction shooting. In this method of shoot- 
ing the compressional waves recorded follow 
essentially the refraction paths characteristic of 
a layered medium. By recording these waves at 
various distances from the shot point and 
plotting suitable time-distance curves the wave 
velocities in high speed layers such as limestone 
may be determined, and also their depth beneath 
the surface. All rock formations may be identi- 
fied by their characteristic velocities of propaga- 
tion of sound waves. 

The reflection record actually obtained, even 
with all precautions for selective filtering of 
certain characteristic frequencies, is very com- 
plex because of the large number of lithologic dis- 
continuities, and one of the major problems of 
the interpreter is to identify corresponding 
reflections from records taken some distance 
apart at the earth’s surface. The refraction 
records are also complex. Energy passes from the 
shot point to the detector by many paths, and 
after the first arrival, energy is received by the 
recording system for several seconds with a 
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shot-detector interval of a few miles. While the 
energy arriving by any one path would be re- 
corded as a relatively short wave train, the 
superposition of a multiplicity of wave trains, 
coming by various paths, makes it difficult to 
identify the individual wave trains. Their 
identification however is the problem which the 
interpreter must solve. Needless to say this 
general type of work requires a basic knowledge 
of the theory of elasticity and the theory of 
vibrating or oscillatory systems. 

The use of geophysics for discovering petro- 
leum reserves is itself a major industry developed 
entirely by the physicist. The oil companies are 
spending at the present time more than 1.5 
million dollars per month on this type of indus- 
trial physics. 

So much for the part physics plays in the dis- 
covery of an oil field. We now come to the 
drilling of the first well. The use by the engineer 
of basic principles of physics has accomplished 
much in perfecting the heavy equipment neces- 
sary to drill a well a mile or two in depth. The 
present record is held by the Gulf Production 
Company at approximately 23 miles, but this is 
by no means the economic limit. The drilling of 
a straight hole was a serious problem five years 
ago and many very ingenious physical instru- 
ments have been developed for surveying well 
bores. However, by the use of proper instruments 
for pressure control on the bit, holes are now 
seldom drilled out of plumb by more than 1 or 2° 
at any point. 

Merely as a detailed example of the use of 
colloid physics in dril'ing, a single application 
to drilling mud will be cited. In the rotary drilling 
of a well, mud is circulated down the drill pipe, 
through the bit, and up to the surface between 
the pipe and the wall of the well. The mud 
column carries out the cuttings and drops them 
in a settling pit after which it is recirculated. 
The mud also serves to plaster the hold prevent- 
ing seepage of well fluid into surrounding forma- 
tions. The mud must be a stable colloidal 
medium in which noncolloidal particles may be 
suspended to increase the density, thereby fur- 
nishing a hydrostatic head sufficient to prevent 
flow of fluid into the well during drilling opera- 
tions. Formation pressures of several thousand 
pounds per square inch are encountered which. 
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if not held in check by the drilling mud pressure, 
may effect disastrous results. The colloidal sus- 
pension is stabilized by negative electrostatic 
charges adsorbed on the particle surfaces. 
Various chemical and physical means are em- 
ployed to secure the proper charge for maximum 
stability. These weighted drilling muds in ex- 
treme cases have a density of 2, and in common 
practice a density from 1.2 to 1.5. In spite of the 
heavy mixture the mud must possess a low 
dynamic viscosity such that it is easily handled 
by pumps without excessive pressure, and such 
that it may readily release cuttings in the settling 
pit. By a low dynamic viscosity is meant the high 


fluidity which is imparted to a thixotropic sus- . 


pension by work inducing shear. These fluids 
are non-Newtonian in that their measured 
viscosity depends decidedly upon their rate of 
shear, the shearing stress to which they are sub- 
jected, and the time during which the shear has 
been effective. As a practical example, if trouble 
is encountered with drilling machinery and the 
well must be shut down, the fluid sets as a gel 
imprisoning the cuttings without settling. Set- 
tling of a mile column of cuttings might result 
in the loss of the well and drilling tools. When 
work is again started, the shearing action of the 
rotary drill stem on the mud immediately breaks 
down the gel into a fluid of low viscosity which 
can be handled by the pumps. 

The next phase in the petroleum industry is 
production, and here the physicist who has 
specialized in hydrodynamics has his oppor- 
tunity. 

It has been necessary to develop a new hydro- 
dynamics of the flow of homogeneous fluids 
through porous media to replace the classical 
Stokes-Navier hydrodynamics. The understand- 
ing of the flow of oil or gas through the under- 
ground reservoir rocks into the producing wells 
was entirely impossible from a practical point of 
view, on the basis of the classical hydrodynamics, 
in which it is assumed that the whole of the 
region of flow is permeated by the fluid. This 
has necessitated the development of essentially 
new fundamental physical laws, in which the 
honeycombed rock structure, composed of im- 
pervious particles connected in irregular net- 
works, is considered, together with the fluid, as 
an integral part of the flow system. When these 
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new dynamical laws are supplemented by the 
physical principle of the conservation of matter 
and by the thermodynamic characterization of 
the nature of the fluid, a self-contained and 
complete new system of hydrodynamics is ob- 
tained which provides the physical mechanism 
for the understanding of the nature of the 
homogeneous fluid motion through porots media. 

The development of a hydrodynamics of the 
flow of heterogeneous fluids—gas-liquid mixtures 
—through porous media is still in its infancy. 
The classical physics gave no consideration to 
the simultaneous flow of gases and liquids even 
through vessels free of sand or other porous 
media. Yet the fundamental physical problem 
associated with the production of oil is that of 
the flow of such mixtures through the tortuous 
and irregular pores of the reservoir rocks. While 
the study of the hydrodynamics of homogeneous 
fluids through such porous media represents a 
natural foundation for the investigation of 
gas-liquid mixture flow, the latter actually repre- 
sents a radically different type of physical 
problem which can be treated only by physical 
methods of analysis. 

As a practical result of this academic work we 
are beginning to understand the mechanics of 
the physical production of oil. Such specific 
problems as water coning, water encroachment, 
water flooding, acid treatment of lime forma- 
tions, may be attacked in an intelligent manner. 

Instruments have been developed for securing 
a sample of well fluids, gas and oil, at the bottom 
of a deep well under actual formation pressure, 
and simple means have been devised for analyz- 
ing the nature of these fluids when brought to the 
surface. Extremely accurate pressure gauges 
have been constructed which permit the record- 
ing of pressures at any depth in a producing 
well. Aided by data of this character the physi- 
cist is able to interpret the pressure variations 
of a producing oil reservoir and predict its 
future performance. Recommendations can be 
made for the most efficient method of production 
so as to obtain maximum recovery of the oil. 
These data also permit estimates of the total 
oil present in a reservoir, and furnish the only 
means for a scientific economic valuation of a 
property. 

By understanding the hydrodynamics of the 
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flow of oil-gas mixtures through tubing, it is 
possible to prolong the flowing life of a well and 
thus delay the installation of costly pumping 
equipment. A large portion of the oil produced is 
in the form of a water-in-oil emulsion. An im- 
portant production problem is the breaking of 
these emulsions by electrical or chemical means. 
There is plenty of physics involved even in the 
latter method, and colloidal physics or surface 
physics, usually referred to as surface chemistry, 
is a fundamental adjunct to petroleum physics. 

We must be content with a single further 
example which is one of our newest developments. 
On the basis of our knowledge of the hydro- 
dynamics of the flow of oil, gas and water into a 
well, it has been possible for our chemists to 
effect a simple method for shutting off water 
horizons without damaging the formations con- 
taining oil. Silicon tetrachloride is injected under 
pressure for a distance of a few feet into the 
producing formations. This reacts with the 
water, forming a precipitate which decreases the 
permeability of the water sand to a very low 
value. Practically 100 percent water shut-off 
has been obtained in preliminary field experi- 
ments. The method may have important com- 
mercial possibilities. 

Petroleum and its products are transported by 
pipe line, by tank ship, by tank car and by tank 
truck. In all of these divisions of transportation 
the physicist plays a role. An obvious applica- 
tion of physics is in the flow of fluid through 
pipe lines where conditions may be either 
viscous or turbulent. Correlation of empirical 
data on this subject by applying dimensional 
methods and the Reynolds’ criterion has greatly 
improved accuracy in design and efficiency of 
operation. A specific field is the investigation of 
surge conditions over a long line, especially in 
starting, stopping, or changing the rate of flow. 
Even the instruments used for controlling opera- 
tions such as flow meters, pressure gauges, and 
tank level indicators owe most of their improve- 
ment in accuracy to physics. 

Refining is generally considered as a series of 
chemical processes, but successful operation is 
largely dependent upon the correct application 
of fundamental physical principles. Not only does 
physics provide a sound foundation on which the 
technology of chemical refining is based, but in 
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many refining operations physical rather than 
chemical changes are involved. 

Distillation, the oldest as well as the most 
important refining process, is, in its simplest 
case, of a purely physical nature. The principles 
involved have long been familiar to the physicist. 
The advent of pyrolytic distillation or ‘“‘crack- 
ing’ with the development of complicated 
modern cracking stills has been made possible 
only through a much wider application of 
physics. The thermodynamics of complex hydro- 
carbon mixtures, liquid-gas phase relations at 
high pressures and temperatures, the viscosity 
of fluids and the conditions for turbulent flow, 
heat transfer by radiation and convection—these 
are some of the problems encountered in cracking 
which the physicist has successfully solved. 

Adsorption is another physical process which 
has been universally adopted as a refining 
method. Removal of undesirable constituents 
from lubricating oil and petroleum distillates by 
means of finely divided clays and other adsorb- 
ents has been shown to follow well-known 
physical laws, the application of which has 
placed these processes on a commercial basis. 
Units for treating over 100,000 barrels of gaso- 
line daily by vapor phase adsorption methods 
have now been installed in this country and 
abroad. 

Solvent extraction of petroleum products is a 
physical refining method which is becoming 
increasingly important. By the use of selective 
solvents less desirable classes of hydrocarbons 
can be eliminated from petroleum raffinates 
making possible the production of premium 
products from ordinary or even inferior crudes. 
A thorough knowledge of distribution coefficients 
and the physical principles of extraction was re- 
quired for the development of these processes. 

Crystallization of wax in readily filterable form 
is an important physical problem in the refining 
of lubricating oil. Application of the results of 
the physics of crystal growth has shown how it 
is possible to avoid ‘“‘shock-chilling’’ and gela- 
tion. Physics also played a considerable part in 
the successful solution of the opposite problem, 
that of preventing wax separation from lubricat- 
ing oils when used at low temperatures. 

A very important contribution of physics for 
which the refining industry as well as all other 
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branches of modern industry is indebted is the 
development of precision instruments for auto- 
matic control. Many of the major oil companies 
are employing such instruments to regulate the 
quality of products obtained in unit processes 
such as distillation or filtration. 

And now for a word about marketing, the last 
phase in the sequential picture we are presenting. 
Two major marketing outlets for petroleum 
products are the service station and direct sales 
to large industrial customers. The choice and 
lay-out of equipment for dispensing gasoline and 
oil are largely based upon the same general 
physical information and methods described for 
pipe lines. In industrial sales the solution of most 
lubrication problems presented by the customer 
is based almost entirely upon the laws of ‘hy- 
draulics, thermodynamics, and the physical 
properties of the lubricant. Marketing with its 
able colleague advertising, is becoming more 
scientific every year. The many problems of 
meeting competitive sales, of breaking down 
consumer resistance and of building up good 
will are being attacked with the same type of 
psychological approach that is employed by the 
experimental physicist. In fact these business 
enterprises are using our ideas, our methods, our 
tools, and equipment, and they are beginning to 
employ our fellow physicists as a necessary 
adjunct to their staffs. 


FOOTE 


And, finally, I must mention the work which 
the physicist conducts in the study and control 
of the physical properties and use of the manu- 
factured petroleum products. The functioning 
of gasoline in an internal combustion motor 
involves such problems as atomization in the 
jet, transportation of droplets in a gas, and 
acoustic properties of the pressure wave in the 
combustion chamber, particularly under condi- 
tions of detonation. Thermodynamics is of 
course the guiding factor in the behavior of the 
combustible mixture in an actual engine. Lubri- 
cation until recent years was almost entirely an 
empiric branch of engineering. In the last decade 
however research in lubrication has been uni- 
versally based on the development of rational 
physical theory for which empirical test data 
serve purely as confirmation. But lubrication 
is only one of the 200 different generic uses of oil. 
Different oils are required for conditioning such 
fibers as cotton, hair, hemp, jute, rayon, silk, 
wool, and other materials such as_ leather, 
feathers, fur, paper and tobacco. Suitable pe- 
troleum products must be devised for new 
processes of absorption, air-cleaning, heat trans- 
fer, flotation, gas fogging and special electrical 
developments. The list of new materials and 
uses is endless, and in a large majority of develop- 
ments the experimental physicist has assumed 
and will continue to assume an important role. 
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New Electronic Tubes and New Uses* 


DAYTON ULREY, Westinghouse Electric and Manufacturing Company, East Pittsburgh, Pennsylvania 
(Received January 23, 1936) 


EVERAL years ago it was probably generally 
believed that the only uses of electronic tubes 
of commercial importance would be, in the future, 
as it had been in the past, in the radio field. 
However, the new applications which have been 
found for so-called ‘‘old type’’ tubes and the new 
types of tubes which have been developed for use 
in entirely new fields during the last few years 
make it appear that this view requires serious 
modification. 


X-Ray TuBEsS 


As illustrative of the importance of a new use 
for an old type of tube we may consider the 
recent developments with x-rays. The x-ray 
tube is one of the oldest of electronic tubes. Its 
value in the past in the medical field for radi- 
ographic, radioscopic, and therapeutic purposes 
has been inestimable and this value will ap- 
parently further increase. Then there are the 
more recently developed and developing applica- 
tions in industry for the examination and analysis 
of metals.! But there is also another field of great 
scientific, and possibly later commercial, impor- 
tance. The development of new species of plants 
and animals, not long ago a relatively rare oc- 
currence and totally beyond the influence of man, 
is now being induced by x-radiation. Mutations of 
wide variety, that is variations which are repro- 
duced in succeeding generations, have been pro- 
duced in many plants and animals. For example, 
over 800 different mutations have been induced 
in barley alone by x-radiation of the seed, and 
even greater numbers of variations of Drosophila 
(fruit fly) have been produced by x-radiation of 
one of the parents. 

A type of variation produced in barley is 
illustrated in Fig. 1 where the plant on the left 
shows the normal growth while the vine-like one 
on the right is a mutation grown from x-irradiated 
seed. 

*From a paper presented at the Conference on Applied 
Physics held at Pittsburgh on November 15, 1935. 

*George L. Clark, ‘‘X-Rays—What Should We Know 


About Them,” Elec. Eng. 54, January (1935). 
+ Taken from Sinnett and Dunn, Principles of Genetics. 
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Fic. 1. Mutations of barley. 


As an incidental—but an important one—in 
mutation experiments with x-irradiated seed is 
the abnormal rate of growth found with some 
plants from such seed. For example, gladioli 
grown from irradiated seed have shown more 
development in two years of growth than is 
found under otherwise similar conditions in three 
to four years of normal growth. Professional 
breeders of plants and animals have thus the 
means of developing more new species in a few 
short years than was formerly possible to obtain 
in a lifetime. One’s imagination alone limits one 
in forecasting the ultimate importance of this 
new field of application for x-ray tubes. 


CATHODE-RAY TUBES 


The cathode-ray tube, though old in principle, 
has been brought to life by improvements of the 
last few years. These have made possible the 
practical use of such tubes for all kinds of oscillo- 
scopic work and no laboratory where circuit 
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Fic, 2. Structure of 6E5 cathode-ray tube. 
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Fic. 3. Typical circuit and appearance of screen of 6E5 
cathode-ray tube. 


studies are made is now considered complete 
without at least one such instrument. Advances 
in fluorescent material make it possible to 
provide a wide variety of color in the light from 
the screen so that the different needs, if not even 
the fancies of different investigators, may be met. 
By the use of screens giving blue and violet light 
it is now possible to photograph transient phe- 
nomena whereas formerly only those phenomena 
which were accurately recurrent could be photo- 
graphed. Other screen materials which have a 
characteristic persistence of luminescences are 
used to prevent ‘‘flickering’’ when slowly re- 
current electrical changes are to be viewed. 

An interesting modification of the cathode-ray 
tube is found in the new electron-ray tube, 6E5, 
in which a visual measure of the voltage im- 
pressed on the grid is given by the angle or area 
of fluorescence on the screen. The structure, 
typical circuit and characteristics are shown in 
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Fic. 4. Characteristics of cathode-ray tube 6ES. 


Figs. 2, 3 and 4. This tube has been used, so far 
as the author is aware, only as a tuning indicator 
for radio receiving sets but since a change of one 
volt on the grid makes a change in the shadow 
angle of almost 20° and since the grid takes 
practically no current, the tube seems to lend 
itself admirably for use as a visual indicator in 
“null’’ circuits. 

Of course, the principal new application for 
cathode-ray tubes is in television. This latent 
application has obviously been delayed by the 
depression but with returning prosperity and 
with active steps being taken in at least four 
different countries to begin commercialization of 
this project, we may expect soon to see a sudden 
growth in demand for such tubes. 


PHOTO-TUBES 


The application of photo-tubes in industry has 
developed very rapidly in recent years because of 
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a combination of circumstances. First, the sensi- 
tivity has been. increased by about 100 to 1; 
second, the cost has been reduced by almost 10 
to 1; and, third, suitably dependable, relatively 
low priced, high vacuum amplifier tubes have 
been made available for use with them. This has 
encouraged literally hundreds of different appli- 
cations? such as counting numbers of automatic- 
machine-made units, automatic traffic control, 
sorting of beans, coffee, cigars, etc., door openers, 
temperature control, elevator control, paper 
cutter control, color matching, transparency 
metering, concentration control, illumination 
control, burglar alarms, and many others. 
Ingenuity has still not been exhausted in this 
direction as evidenced by the fact that every 
issue of Flectronics describes still others worthy 
of mention. Practically all of these are, however, 
essentially similar in principle, consisting of a 
light source, photo-tube with associated ampli- 
fier, and suitable relays. 

Quite a different type of application has lately 
been found. The modern demands for artificial 
sources of ultraviolet light for radiation therapy, 
for bacteria sterilization, and for irradiation of 
foods to improve the vitamin content have also 
created a demand for means of accurately meas- 
uring this ultraviolet dosage. This necessitates a 
photo-tube insensitive to radiation outside the 
particular range of interest. For example, radia- 
tion most valuable for therapeutic purposes is 
believed to consist only of the narrow band 
between 2800 and 3100 Angstroms. To measure 
this without interference from longer or shorter 
wave-lengths in the light source, a photo-tube is 
used having a titanium cathode, the long wave- 
length limit (the longest wave-length for which 
it is photoelectrically active) of which is 3200A. 
To cut off the other end of the ultraviolet a type 
of glass is used for the tube which is entirely 
opaque for wave-lengths below 2700A. Thus the 
tube is responsive only to wave-lengths between 
2700 and 3200A. Similarly, radiation most useful 
for bactericidal purposes is considered to lie 
below 2800A. A photo-tube for measurement of 
this radiation employs a tungsten cathode, the 
long wave-length limit of which is 2700A, and a 
glass is used which absorbs all radiation below 


?See V. K. Zworykin and E. D. Wilson, Photocel!s and 
Their Applications, second edition (John Wiley & Sons). 


2200A. The tube, therefore, is sensitive only to 
radiation between 2200 and 2700A. 


Low Grip CURRENT TRIODES 


For the proper use of these photo-tubes, as 
well as for other applications where minute 
currents or voltages are to be amplified, special 
high vacuum amplifiers of unusually high grid 
resistance are required. Recently a line of such 
amplifier tubes has become available so that the 
engineer can now design for such applications 
with due consideration to both cost and per- 
formance. Typical tubes of this type are shown in 
Fig. 5 while in the table below are indicated the 
typical order of grid impedances and types of 
applications for each. 


HIGH FREQUENCY OSCILLATORS 


Other new types of high vacuum triodes, not 
differing from those just mentioned in basic 
principles of structure, but radically—even spec- 
tacularly—different in shape, in size, and in type 
of application, are those developed for the genera- 
tion of so-called “‘ultra-high”’ frequencies, that is, 
frequencies of the order of one hundred times 
those used in radio broadcasting. 

One of the most interesting and important 
applications of such tubes is in the field of 
physical therapy. High frequency oscillators are 
used for this purpose, in general in two different 
ways. First, equipments developing an output of 
one to ten kilowatts are used to develop a general 
fever with temperatures of 104 to 106°. This 
treatment has been definitely proven to be 
effective in combating such diseases as arthritis, 
bursitis, multiple sclerosis, paresis, and others. 
Other means of producing artificial fevers are now 
available, however, and apparently give equally 
good results and consequently this application 
does not assume the major importance it other- 
wise would have. 

The other type of use does not involve appreci- 
able elevation of the body temperature but an 
output of 150 to 250 watts at about 50 million 
cycles is applied locally to a relatively small 
affected area. An enormous amount of investiga- 
tion, particularly by Schliephake* in Germany, 


3E. Schliephake, Short Wave Therapy (The Actinic 
Press). 
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Fic. 5. 


Type of tube Grid impedance 


Typical amplifier use 


Radio receiving tube 10° ohms (plate V = 100) Radio and audiofrequencies. 
RJ-553 10® ohms (plate V = 100) Weak photo or similar current. 


RH-506 10'? ohms (plate V =6) 
FP-54 10'6 ohms (plate V =6) 


in the treatment by this method of boils, car- 
buncles, bone tumors, focal infections, pulmonary 
abscesses, etc., has shown almost miraculous 
results. For example, a survey of international 
literature indicates that the mortality in the case 
of pulmonary abscesses has been 60 to 75 percent 
with conservative treatment and 30 to 40 percent 
for surgical treatment. Schliephake has, however, 
treated 20 different cases, not specially selected, 
without a single death. 

It should be added that Schliephake claims to 
have proven conclusively that in the treatment of 
most of such disorders, at least, there is a definite 
selective effect. For example, in one of his papers 
he says “It is a striking fact in this connection 
that wave-lengths of four to six meters worked 
with the certainty of an experiment, while wave- 
lengths of eight to fifteen meters, in spite of a 
greater output, were ineffective.”’ 

Results such as these have stimulated great 
activity by the medical profession in this country 
and it is not improbable that this type of equip- 
ment will soon be considered fully as essential in 
hospitals and physicians’ offices as x-ray ap- 
paratus now is. 

Equipment similar to that last described is also 
being used by biologists in studies on plants and 


Extremely minute photo-currents, H2 ion concentration, etc. 
Shot effect, etc. 


animals and some evidence, though possibly no 
proof as yet, has been obtained that this treat- 
ment also induces mutations. 

There are probably numerous other applica- 
tions for power at these high frequencies when 
this is available commercially in sufficient magni- 
tude and when sufficient experimental work has 
been done to establish values. As an example, 
may be mentioned the treatment of grains and 
cereals for insect extermination. Rather extensive 
experiments in the Westinghouse Laboratories 
and by the B & O Railroad have shown that 
weevil infested wheat can be completely de- 
infested by this means at a cost of not more than 
a few cents per bushel. In these experiments, 
wheat was passed between condenser plates in 
the oscillating circuit of an oscillator delivering 
20 to 30 kw at five to six meters wave-length. The 
adult insects, their eggs and larvae each have 
relatively high dielectric losses compared to the 
wheat and consequently rise in temperature ver 
rapidly—so rapidly in fact that one second ex-| 
posure in the more powerful high frequency field 
proved to be a lethal dose. 

An idea of the size and nature of the weevils 
thus treated can be had from Fig. 6 in which a 


grain of wheat is shown for comparison. 
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Fic. 6. Showing relative sizes of grain of wheat and weevils. 
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Fic. 7. Circuit of standing wave oscillator. 


The circuit used in these experiments, known 
as “the standing wave oscillator,’ is illustrated 
in Fig. 7, the subject to be treated being placed in 
or passed through the treating condenser. The 
tube used is shown in Fig. 8. This combination 
provides an output of 13 or more kw in the range 
of 40 to 70 megacycles. By the use of a second 
AW-200 tube in the other end of the oscillator 
working in “push-pull” with the one shown the 

*1. E. Mouromtseff and H. V. Noble, ‘‘A New Type of 


Ultra-Short Wave Oscillator,’”’ Proc. Inst. Radio Eng. 20, 
August (1932). 


Fic. 8. AW-200 Ultra-high frequency, high power 
oscillator-amplifier. 


output has been approximately doubled, the 
highest power yet developed anywhere at such 
frequencies. 

A view of part of the experimental set-up is 
shown in Fig. 9, this also illustrating some of the 
spectacular phenomena producible with high 
power at high frequency. The operator has 
started the “ionic flame’’ merely by touching 
with a conductor on the end of an insulating pole, 
one end of a rod of half wave-length placed 
parallel to but out of electrical contact with the 
standing wave tube. 

It may be mentioned that essentially this same 
equipment has been used in substantially the 
same manner for other applications such as 
heating—even melting—of glass and other in- 
sulating materials, the rapid roasting of meat, 
etc., where there is an obvious advantage of 
developing heat in all parts of the object simul- 
taneously, rather than having to delay for con- 
duction from the outside surfaces to the center. 
It is possible that applications of this type may 
be commonplace in the future. 

Other applications for tubes developing high 
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_ Fic. 9. Showing experimental set-up of “‘ultra-high” 
frequency high power oscillator and the production of 
‘ionic flame.” 


frequency currents, though not so high as in the 
cases just mentioned, are found in the production 
of high temperatures. For example, the melting 
of the more precious metals, such as gold and 
silver, to form castings or alloys can best be 
done in a high frequency furnace where freedom 
from contamination can easily be obtained by 
providing an inert atmosphere or even a high 
vacuum, if desirable. Such equipment is par- 
ticularly useful in the research laboratory because 
of its great flexibility and adaptability to differ- 
ent problems.’ In the Westinghouse Research 
Laboratories, for example, there are no less than 
eight of these units capable of supplying outputs 
of from 10 to 30 kw at frequencies of from 100,000 
up to 2,000,000 cycles. 


5H. V. Noble, ‘‘Industrial High Frequency Power,” 
Electronics, Oct. (1935). 
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SECONDARY EMIssION MULTIPLIERS 


A new type of electronic amplifier, departing 
radically in structure and mode of operation from 
those mentioned before, is now being developed 
in several laboratories. Within the last year 
Farnsworth described, before the Pittsburgh 
Physical Society, a type of amplifier which has 
been termed a “‘secondary-emission multiplier’ 
and within the last month Zworykin has de- 
scribed before the Institute of Radio Engineers in 
New York a quite different variety of tube,® also 
developed on the secondary-emission principle. 

The principle of this type of amplifier is 
illustrated in Fig. 10. Electrons are supplied from 
a primary cathode by some means, here indicated 
as photoelectric, are attracted to anodes of 
progressively higher potential by means of com- 
bined magnetic and electrostatic fields. The 
dotted lines indicate the direction but not the 
relative number of electrons for at each anode 
their number is increased by a factor which may 
be as high as six or seven. Thus for a ten-anode 
unit there should be of the order of 100 million 
electrons arriving at the final anode for each pri- 
mary electron from the cathode. Experimentally 
it has been found that this multiplier amplifier 
has the following inherent advantages: 


1. Current amplifications of several million can be obtained 
in a single tube. 

2. Amplification does not vary with frequency as is 
ordinarily the case in multi-stage conventional type 
amplifiers, 

. The signal to noise ratio (that is, noise developed due to 
shot effect in the vacuum tube) is from 50 to 100 times 
better than in conventional amplifiers of equal gain. 


It seems probable therefore that this new type 
of tube may before long find its place in our high 
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Fic. 10. Scheme of an electron multiplier tube. Magnetic 
field is applied perpendicularly to the plane of the drawing. 


6 V. K. Zworykin, ‘‘Secondary Emission Electron Multi- 
pliers,” Electronics, Nov. (1935). To be published in more 
detail in a forthcoming issue of Proc. Inst. Radio Eng. 
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gain amplifiers such as are found in radio broad- 
cast transmitters and as will be required par- 
ticularly in television transmitters—not im- 
probably also in sound and television receiving 
sets. 


Gas TRIODES 


The tubes thus far mentioned are all of the 
high vacuum type. While the use of high vacuum 
provides characteristics which are essential to the 
operation of such devices for the purposes for 
which they are used, it also involves limitations 
which practically preclude the possibility of their 
use for many other purposes. These principal 
limitations are: (1) A relatively low current 
capacity; (2) the necessity for relatively high 
anode voltages; and (3) the high internal losses 
resulting in low over-all efficiencies. All of these 
undesirable characteristics, particularly in tubes 
used for power and control purposes, can be 
largely eliminated by the use of gas or vapors at 
low pressures instead of a high vacuum but at the 
cost of losing the property of linear amplification, 
or, in other words, of having the anode current 
proportional to the control grid voltage. The 
_ grid in a gas triode is effective only in initiating a 
current between the cathode and anode but 
after the discharge has once started further in- 
crease or decrease of the grid voltage has no 
effect whatever on the anode current. The cur- 
rent can, therefore, only be interrupted by break- 
ing the circuit or reducing the anode voltage to a 
value too low to sustain the discharge. From such 
characteristics it will be seen that these gas tubes 
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Fic. 11. Showing structure of a gas triode. 
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Fic. 12. Characteristics of the gas triodes. 


lend themselves well for use in a.c. circuits where 
the discharge must be restarted every cycle. The 
phase or magnitude of the grid voltage will then 
determine at what part of the cycle the discharge 
is started and therefore the fraction of the total 
half-wave that passes through the tube each 
cycle. A typical structure of this type of tube and 
typical characteristics are shown in Figs. 11 and 
12, respectively. Time does not permit of de- 
scription in detail of any of these tubes, known 
variously as gas triodes, controlled rectifiers, 
thyratrons, grid glow tubes, trianodes, etc., nor 
of detailed manner in which they are applied in 
so many different ways in industry.’ Only a few 
examples of these applications may be men- 
tioned: There are innumerable automatic ma- 
chine operations where the gas triode acts as a 
high speed, noiseless, contactless, reliable relay 
requiring no power to excite it; they are used in 
the control of theater lighting making it possible 
not only to displace many bulky, power con- 
suming rheostats, but also to pre-set as many as 
ten different lighting arrangements, many of 


7 See Gulliksen and Vedder, Industrial Electronics (John 
Wiley & Sons); and Henney, Electron Tube in Industry 
(McGraw-Hill Publishing Co.). 
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which were previously entirely impracticable, 
and still be able to fade smoothly from one scene 
to another as the show progresses; they are used 
for the automatic regulation of voltage from 
large generators; they provide a means of con- 


verting d.c.power to a.c.power and vice versa at. 


high voltages thus making possible the d.c. 
transmission of power; they are used for the 
speed control of motors for high speed photog- 
raphy, for telemetering, etc. These and many 
others indicate the trend of industry to adjust 
itself rapidly to the new possibilities of this 
relatively new tool. 

The thermionic gas triode, to which the fore- 
going description applies, is available in sizes 
capable of delivering currents of from the order 
of a hundred amperes down to small values. 
There are, however, many applications which 
require enormously high currents for short 
periods even though the average current required 
is small. For such use the thermionic cathode 
must be large enough to supply the peak demand. 
This is impractical where the peak currents are in 
the order of thousands of amperes and in such 
cases pool type cathodes with keep-alive arcs 
have been resorted to though at the expense of 
considerable complication. 


IGNITRONS 


A recent development of a new type of tube,*® 
which has been called the ‘‘Ignitron’’ because of 
the manner in which the discharge is started, has 
overcome this former limitation of low ratio of 
peak to average current and opened still further 
the door to the application of tubes to industry. 

The principles of structure and operation of 
this device are shown diagrammatically in Fig. 
13. A rod of suitable resistance material, such as 
boron or silicon carbide, is partly immersed in the 
mercury pool. This has the property that when a 
suitable current is passed through this mercury- 
resistor junction a minute spark is formed at the 
surface of contact and this acts as an initiating 
cathode spot permitting a discharge to pass if 
anode voltage is applied. In this type of circuit 
the point of breakdown, e,, is determined only by 


’ J. Slepian and L. R. Ludwig, ““A New Method of 
Initiating the Cathode of an Arc,” Trans. A. I. E. E. 52, 
(1933). D. D. Knowles and E. G. Bangratz, ‘‘The Igni- 
tron,” Electric J., December (1933). 
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Fic. 13. Structure and operation characteristics of the 
ignitron. 


the time required for the anode to reach break- 
down voltage, assuming that the resistance R is 
low enough to allow sufficient current to flow 
through the ignitor at this voltage, and anode 
current then flows for the remainder of the half- 
cycle. The average anode current may however be 
controlled, as in a gas triode, by varying the time 
of current starting. The method of accomplishing 
this is illustrated in Fig. 14. Here the anode 
voltage rises above the critical breakdown value 
but current does not start until the current in the 
ignitor circuit reaches a critical value. This in 
turn is controlled by the bias on the grid of the 
gas triode. It is obvious that by this means any 
predetermined portion of the positive half-wave 
may be passed each cycle. 

Another type of firing circuit particularly 
suited for precision work is illustrated in Fig. 15. 

As an indication of the type of ratings which 
ignitrons may have, Fig. 16 shows one built in a 
ten-inch glass blank which can safely pass surge 
currents of 2000 amp. though the average current 
rating is 25-50 amp. 

Other sizes, up to average current ratings of 
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several hundred amperes in all-steel structures, 
are also available. 

Such tubes, though less than three years old 
have already found extensive uses in industry, 
particularly in spot and seam welding, such as is 


USES 


Fic. 16. Ignitron. Average current rating 25 to 50 amperes. 
Surge currents 2000 amperes. 


used in fabrication of “‘all metal” tubes, of light 
weight, high strength structures for bridges, 
airplanes, railway cars, etc. Many other applica- 
tions are in sight where accurately controlled high 
current surges are necessary. | 

The foregoing discussion, though obviously 
quite superficial and incomplete, is perhaps 
sufficient to indicate the wide variety of relatively 
new applications of electron tubes in every day 
use. This is in spite of the fact that most of this 
development has grown in depression years. 

With better business conditions and with the 
gradual absorption by industry of college gradu- 
ates who have been given instruction in the 
fundamentals of electron tube properties and 
uses, one may forecast even more rapid develop- 
ments in this field in the next few years. 
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The Theory of Acid Treatment of Oil Wells Producing from Limestone Reservoirs 


M. Muskat Anp R. D. Wyckorr, Gulf Reseerch and Development Corporation, Pittsburgh 


(Received November 25, 1935) 


A theory is given of the acid treatment of oil wells 
producing from limestone reservoirs, based upon the as- 
sumption of a potential flow of the liquid in the limestone. 
Calculations are presented for the increase in production 
capacity of a well resulting from an increase of the permea- 
bility in an annular ring surrounding the well bore, obtained 
by the reaction of the acid (HCl) with the limestone. It is 
found that the effects are maximal if this inner zone is 
initially of very low permeability, and is increased to a 
normal value, while if the permeability is normal originally, 
the effect of the acid treatment will be small. For limestones 
that are highly fractured, a representation is used wherein 


the fracture is replaced by a linear porous medium, of 
width equal to the fracture, bisecting the limestone proper 
and of a permeability that is equivalent to its fluid carrying 
capacity as given by the classical hydrodynamics for a free 
linear channel. The effect of the acid is then taken to be 
that of widening the fractures, thus increasing their effec- 
tive permeability. The analysis here shows that increases in 
the production capacities as high as one hundred-fold can 
be obtained in such fractured limestone systems for 
increases in the fracture width no greater than 1 mm, the 
effects being greatest when the initial fracture widths are 
small. 





I. INTRODUCTION 


— treatment” is the name given to the 
process of artificially stimulating the 
production of oil from a well by the introduction 
of acid into the well bore. The method has thus 
far been found successful only in the case of 
wells producing from limestone reservoirs, the 
acid used being hydrochloric acid, which reacts 
with and dissolves the CaCO; of the limestone.! 
The immediate effect of such a ‘‘treatment”’ will 
evidently be a “‘cleaning’’ of the exposed face of 
the producing horizon. However, if the limestone 
is highly fractured, as is frequently the case, the 
acid may not only penetrate and widen the frac- 
tures, but it may open up previously untapped 
parts of the reservoir. In the former case only the 
production capacity of the well, and hence the 
economic ultimate recovery of oil from the 
producing zone will be increased, while in the 
latter the physical ultimate recovery will be 
increased as well. As the occurrence of and mag- 
nitude of the latter effects will depend largely on 
chance, they will be excluded from the present 
discussion. 

Although it has generally been considered that 
the flow of homogeneous fluids through lime- 
stones cannot be treated in the same manner as 
that through sandstones, because of the irregular 
occurrence of fractures in the former, which 


' Hydrofluoric acid has been tried with wells producing 
from sandstones, but with very little success. 


undoubtedly play a large part in the flow charac- 
teristics, it nevertheless appears to be reasonable 
to apply such treatments in two important limit- 
ing cases. These are: (1) that in which the 
fractures are of limited extent and disseminated 
throughout the main body of the limestone with 
some degree of uniformity, i.e., a network of 
small fractures; and (2) that in which the lime- 
stone as a whole is largely free of fractures, except 
for a limited number of extended length. The 
macroscopic features of the former will evidently 
be quite the same as those of an ordinary homo- 
geneous sandstone, with an effective permeability’ 
which will clearly be considerably higher than 
that of the limestone proper. In the latter case, 
the flow system may be considered as a compo- 
site synthesis of a uniform medium representing 
the limestone proper bisected by others of differ- 
ent permeability, corresponding to the fractures. 
While these fractures are, in principle, flow 


‘channels free of any porous medium, they may be 


given “‘effective’’ permeabilities to correspond to 
their carrying capacity, as determined by the 
classical hydrodynamics. As these effective 
permeabilities will be very much higher than that 


2 The permeability of a porous medium is the volume of 
fluid of a unit viscosity passing through a unit cross section 
of the medium in unit time under the action of a unit 
pressure gradient. In the c.g.s. system, with the units of 
pressure and viscosity taken as the atmosphere and centi- 
poise, the permeability unit has been named the ‘‘darcy.” 
For a detailed discussion of the measurement of this quan- 
tity, cf. R. D. Wyckoff, H. G. Botset, M. Muskat, and 
D. W. Reed, Rev. Sci. Inst. 4, 394 (1933). 
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of the limestone proper, which may be of the 
order of 0.01 darcy, the flow will be essentially 
one in which the limestone proper feeds its fluid 
into the highly permeable crevices or fractures, 
and these bring the fluid directly or by complex 
interconnection into the outlet wells. 

This same distinction between the types of 
reservoir may be made the basis of discussing the 
effect of the acid treatment. Thus we will con- 
sider first the effect of the acid if the reservoir is 
essentially homogeneous—in a macroscopic sense 
—on the assumption that the acid increases the 
permeability of an annular ring concentric with 
the well bore, the system as a whole being radial. 
Then we will treat the case where the limestone is 
highly fractured, the fractures bisecting the 
whole of the main body of limestone. 


Il. THE Errect or Actp TREATMENT 
IN RADIAL SYSTEMS 


The diagrammatic representation of a radial 
limestone flow system, which has been treated by 
’ acid, is given in Fig. 1. The annulus between the 
well bore,’ of radius 7,, and the radius fo is 
supposed to have been affected by the acid, while 
the limestone beyond ro remains with its original 
permeability ke. The resultant effect of the acid 
may be conveniently described by the ratio of the 
production capacity of the well after treatment, 
Q, to that before treatment, Qo. These may be 
obtained from the solution to the potential 
problem corresponding to the system of Fig. 1. 
The formulation of this problem, in turn, follows 
from the application of the equation of continuity 
to Darcy’s law relating the macroscopic vector 
fluid velocity 3 to the pressure gradients in a 
porous medium, namely: 


d= —kV p/n, (1) 


k being the permeability of the medium? and yu 
the viscosity of the fluid. For an incompressible 
liquid, therefore, it follows at once that p must 
be a potential function, i.e. : 


V2p=0. (2) 


*It is here supposed that the limestone is a horizontal 
slab of uniform thickness and that the well completely 
penetrates the producing zone, so that the problem is not 
only two-dimensional, but that furthermore gravity will 
have no effect on the flow within the porous medium. 





acid offected 
zone 


Fic. 1. Diagrammatic representation of a radial flow system 
subjected to acid treatment. 


Supposing, now, that the pressure at the bottom 
of the well bore is kept at p,,, and that the value 
at a large distance r, is p,., the “reservoir pres- 
sure” or “external boundary pressure,” the ana- 
lytical problem may be finally stated as: Find 
two potential functions p; and pe, such that: 


Pi=Pw St Meal 470) 
pi=p2; 


By analogy with the electrical problem of a 
cylindrical condenser, whose dielectric is com- 
posed of two layers of different dielectric con- 
stant, we can write down at once the solutions 
for p; and p2 as: 


Pe— Pw 
pi= 


log 70/fwt+ea log r./ro 


Po= pe > 7=TPLe, 


k,0p;/0r=k20p2/dr : at r=ro. (3) 





log £/fwt Pw; 


¥9="=SPo, , 
(4 
al pe— Pw) : 
p2= ; log r/re+pe; 
log ro/fw+a log r/o 





oS, 


where a=k;/ke. The rate of fluid flow into the 
well, per unit thickness of the porous medium, is, 
by virtue of Eq. (1), given by: 


2rkr Op 2rki(pe— Pw) /p (5) 
= —_——_——_- —- = . | 
un Or log r/rota log r./ro 
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ky /k; 


Fic. 2. Increase in the production capacity of a radial flow system due to 
acid treatment, if the initial permeability is everywhere uniform. Q/Qo = (pro- 
duction capacity after treatment)/(production capacity before treatment). 
k, /k;=(permeability of affected zone after treatment)/(permeability of 
affected zone before treatment). 7>=radius of affected zone; well radius= }’; 


external boundary radius = 500’. 


so that the ratio of the flux in the composite 

system to that, Qo, in the homogeneous system of 
permeability k2(k;=ke, a=1), is: 

Q a log 7, ry ky 

; a=—. 

ke 





(6) 


mest 
Qo log ro ruta log 7, ro 


The difference of Q/Q» from the value unity 
evidently gives the relative change in the pro- 
duction capacity of the well caused by the pres- 
ence of the annulus of permeability k; in the 
immediate neighborhood of the well, and extend- 
ing to the distance rp in a porous medium which 
elsewhere has the permeability ke. 

To show most clearly the implications of Eq. 
(6) with respect to the original problem of acid 
treatment, it is convenient to consider separately 
the cases in which the initial permeability of the 
annular region is normal and equal to that of the 
main body of limestone, so that the acid simply 
increases this permeability to a higher value, 
and that in which the annular region has an 


initial low permeability due to a natural inhomo- . 


geneity in the producing zone or because of 
mudding and plugging in the course of drilling, 
the effect of the acid here being essentially that 
of raising the permeability approximately to its 


normal value. In the first case the effect of the © 
acid treatment is given graphically in Fig. 2, 
where the abscissae k,/k; give the ratio of the 
final annular permeability to its initial value, k; 
being also equal to ke, that of the main body of 
limestone. It will be seen that the effects as a 
whole are not very large unless the radius of the 
affected annulus is quite large, which would re- 
quire very great quantities of acid. Furthermore, 
Q/Qo does not depend markedly upon the abso- 
lute value of the permeability increase as long as 
k,/k; is of the order of 5 or greater. In fact, even if 
the inner annulus is completely dissolved away, 
so that k,/k;= ~, the increase in Q) would equal 
only that corresponding to an increased well bore 
radius from 7,, to ro. Thus the production rate 
will not be doubled unless the effect of the acid 
extends to the radius 11.2’, even though the acid 
completely removes the limestone within that 
radius, for r,=1/4’ and r.=500’. 

The corresponding curves for the second case 
are given in Fig. 3. Here the abscissae are the 
ratios of the initial permeability of the inner ring 
to the value after acid treatment, it being as- 
sumed that this latter value is equal to that of the 
main body of limestone. Large values of Q/Qp are 
here more readily obtained, even for small ro, if 





ACID TREATMENT OF 


WELLS 


-—$—$—$ —_—_>—_—_+—_ 


ki/k, 


Fic. 3. Increase in production capacity of a radial flow system due to acid treat- 
ment, if the acid affected zone is initially of lower permeability than rest of limestone 
and is raised to the normal value by the acid. Q/Qo=(production capacity after 
treatment) /(production capacity before treatment). &;/k: = (permeability of affected 
zone before treatment)/(permeability after treatment). ro=radius of affected zone; 
well radius= }’; external boundary radius = 500’. 


the inner ring is initially of very low permeability, 
and are not limited by the value of 7. As the 
usual treatment, in which some 500 gallons of 18° 
Baume HCI are used in each ‘“‘dose,’’ will pene- 
trate a 10’ thickness of 20 percent porosity 
limestone to a radius of only 3.3 ft., it is clear that 
acid treatment will be of but limited value unless 
the face of the well bore is actually plugged quite 
seriously. For evidently in the limiting case when 
the well face has been completely plugged or 
mudded off, k:=0, and Q/Q) would become 
infinitely large.‘ 


Il]. THe Errect or Acitp TREATMENT IN 
HIGHLY FRACTURED LIMESTONES 


That extended fractures in a limestone reser- 
voir may really play a significant role in the 
production for such reservoirs, is a consequence 
of the fact that a fracture of even a small width 


‘ These results are also of interest in showing the effect 
on the production capacities of wells due to small regions 
of natural inhomogeneity localized about the well bore. 
Indeed, they provide the explanation for the large varia- 
tions frequently observed in the production capacities of 
neighboring and apparently identical wells drilled into 
the same producing horizon. 


may have an effective permeability hundreds of 
times greater than that of the limestone proper. 
For a real fracture of width w may evidently be 
considered as equivalent to an open linear 
channel of equal width. Now the carrying ca- 
pacity of such a linear channel per unit pressure 
gradient for viscous flow conditions is given by 
the classical hydrodynamics by : 


Q=w*/12u, (7) 


where uy is the viscosity of the liquid. The equiva- 
lent permeability of the channel is, therefore: 


=w?/12=108w?/12 darcys, (8) 


if w is expressed in cm. Hence a fracture of only 
0.1 mm width will have a permeability of 833 
darcys, whereas the permeability of the limestone 
proper will usually be of the order of 0.01 darcy. 
In fact, the total fluid carrying capacity of a 
complete radial system of radius 45’ consisting of 
a limestone of permeability 0.01 darcy, can be 
carried by a single linear fracture 45’ long, of 
depth equal to that of the radial system and of 
width 0.126 mm. 

We shall, therefore, base the analytical treat- 
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ment of such limestone-fracture systems upon a 
representation in which the limestone and frac- 
ture are two general contiguous porous media of 
different permeability, that of the fracture being 
the “effective” value given by Eq. (8). Further- 
more, the simplifications will be introduced that 
there is but a single fracture in the system, which 
bisects the otherwise infinitely extended and 
uniform limestone proper, the fracture extending 
throughout the whole thickness of the limestone. 
The well draining the system will be supposed to 
be exposed only to the fracture,® and to penetrate 
the whole depth of the limestone, so as to make 
the problem essentially two dimensional, a sec- 
tion of the typical horizontal plane being shown 
in Fig. 4. 

Although in the practical realization of such a 
system the fracture will be of very small width, 
this fact will be anticipated only to the extent 
that the boundary condition to be imposed at the 
outlet terminal of the medium, the fracture, will 
correspond to the more convenient requirement 
of a uniform flux density which, in view of the 
small fracture width, is quite equivalent to that 
of a uniform pressure, which would be required 
physically. It will also be convenient to take the 
fracture width as the unit of length and both the 
media representing the fracture and limestone to 
be of infinite extent, the latter also being of 
infinite width. Because of the symmetry of the 
system—containing but a single fracture—only a 
single quadrant of the (x, y) plane, the horizontal 
section of the limestone fracture system, (Fig. 4), 
need be considered. The analytical problem may 
then be formulated as: Find pressure distribu- 
tions p; and 2 such that: 


V7p1=0; V*p2.=0, 


dp;/dyv=0 : y=0; dpi/dx=1 : x=0, 


Pi=pe2, kiOpi/Oy=k2dp2/dy : y 


Op2/dx=0, x=0, y>}. 


The resultant values of p; and p2 may be most 
readily obtained by a synthesis of those due to 
unit strength sinks distributed over the y axis for 
0=y=1/2, ie., along OA. For these individual 


5 Because of the very low permeability of the limestone, 


as compared to that of the fracture, the flow into the well. 


directly from the limestone proper may be neglected. 
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solutions, however, we may use those derived 
by one of the authors® for the analogous elec- 
trical problem which is of interest in the theory of 
electrical prospecting. These are: 


a 


© da 
pi=log it (y— hy} +2 f —[(2/6—1)e-« 
0 


+2} (6 sinh (}—h)a+cosh (3—h)a)e-*/A 
+(1—6) cosh hae~*- /A} cos xa], 


© dof l1 
p=4f —|--~ 
0 a 6 


e~-D« cosh y cosh ha cos xa 
1 os 


(10) 





sinh (a/2+7) 


where 6=k2/k,=tanh y; 


(12) 
= —2(5 cosh a/2+sinh a/2), 


it being explicitly assumed in the transformation 
from 6 to y, that 6<1. 

Integrating these expressions with respect to h, 
from h=0 to 1/2, and denoting the integrated 
values also by p: and pe, it is found that: 


bi=(2—y) log {x°+(3—y)*} +y log (x? +’) 


x 
—1+2x tan—' ——_—_——_ 
2x°+2y?—y 


dof /1 1 2 cos xa 
af [Cee 
0 Qa 6 2 aA 
X fe-¥*(6 cosh a/2+sinh a/2—6) 


+(1—6)e-*/?+" sinh a, 2 | 


© da 
= f =| ctnh ye? —2e- aa 
0 a 


cosh y sinh a/2 cos xa 
x 7 | (14) 
a sinh (a/2+y) 





For purposes of computation it is convenient to 
remove the constants by subtracting the values 
of p at x=0. Thus for y=1/2: 


6 M. Muskat, Physics 6, 14 (1935). 
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Fic. 4. Diagrammatic representation of a limestone 
bisected by a fracture of effective permeability ki, with 
the appropriate boundary conditions. 


pi(x, 3) —p1(0, 3) 
=} log (1-+4x*) +2x tan- (1/2x) 
© sin? xa/2da 
vf 
J) o sinh (a/2+7) 
x [e-/* sinh (a/2+7)+sinh (a/2 —v¥)] 








© sin? xa sinh (a—y) 
= xt f da, 
0 a sinh (a+y) 


p2(x, >) — p(0, 2) 


© sin? xa/2 sinh a/2da 
=8 cosh vf ; 
0 a sinh (a/2+y7) 





Several pressure distributions (along y=1/2) 
as given by the above formulas’ are plotted in 
Fig. 5. The pressure variation is linear at first and 
then assumes a logarithmic form, the pressure for 
very large x, and y=10~, approaching the value: 


p=(2 cosh y/y)[1.2704+log yx]. (17) 


However, the linear variation of p persists for 
larger and larger values of x as y is made smaller, 
and in fact becomes strictly linear throughout for 
y=0, when the medium (2) has a vanishing 
permeability. For then it follows from Eq. (15) 
that: 


pi(y =0) =2 rx. (18) 


On the other hand, in the limiting case when 


the medium (2) has the same permeability as 
medium (1), y=, and 


7 For y=10- the integrals of Eqs. (15) or (16) may be 
closely approximated by combinations of simpler ones, 
which may be either evaluated exactly or derived from 
tables of the ‘‘sine’’ and ‘‘cosine integrals.” 


x 1o* 


Fic. 5. Pressure distribution along a fracture fed by a limestone reservoir. 
ki/k2= (effective permeability of fracture)/(permeability of limestone proper). 
Ap=pressure drop from well terminal of fracture to distance x (measured in units of 
fracture width) along fracture, for a production rate of 4rk;/u units per unit thickness 


of limestone. 
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Fic. 6. Resistivity of fracture-limestone systems for unit 
effective permeability of the fracture (&,). Curve I: pressure 
drop, Ap, between the well and a distance 4-10° (in units 
of fracture width) from the well, for a unit flux per unit 
limestone thickness. Ap is in atmospheres when 2; is 1 
darcy. Curve II: fraction of total flow, Q;/Q, entering the 
well which has traveled through the fracture for the whole 
distance of 4-105 units. k2=permeability of limestone 
proper; liquid viscosity = 1 centipoise. 


p(y= ~) =log (1+x7)+2x tan (1/x). (19) 


Observing that Eq. (17) gives the pressure 
differential between the well and the points 
along the fracture at the distance x for a flux of 
2xk,/u into the well from either side of the y axis, 
the pressure differential for unit production rate, 
or the effective resistance of the system, for the 
range of k2/k, of practical interest— > 10°—, and 
for x=4-10°, is plotted as curve I in Fig. 6 for 
k,/u=1. The decreasing slope for decreasing 
ko/k, in Fig. 6 is clearly due to the fact that as the 
permeability of the second medium (limestone) 
continues to decrease, more and more of the 
fluid is carried directly by medium (1)—a fracture 
in the practical realization of the system under 
consideration—so that the effective conductivity 
of the combined system will not be reduced by as 
large a factor as will that of medium (2). In fact, 
it follows from Eq. (16) that: 

A a 
a 4 cosh ¥ 
Ox 


sin xa sinh a/2da 


a sinh (a 2+) 


ctnh ¥ 1 ee 
= | 2-—(enk? »-1) +0(—) | (20) 


x x x* 
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Assuming then that 0p/dx is uniform for large x 
over the width of the fracture, it gives the flux 
flowing directly along the fracture at the large 
values of x. Dividing by 27, one obtains the 
fraction of the total flow into the well which has 
come directly through the fracture, at least be- 
yond the distance x from the well. The curve 
showing how this fraction varies with k,/k:s for 
x =4-10%, is given as curve II in Fig. 6. Thus it is 
seen that whereas for k,/k2=10*, 99.2 percent of 
the total flux is fed into the fracture from the 
limestone within the distance x, only 54 percent 
of the total flux is so derived if k;/k2=10°. 

As curve I of Fig. 6 is based on a fixed (unit) 
permeability of the medium 1, it gives the varia- 
tion of the composite resistance of the limestone- 
fracture system in which the fracture is kept 
fixed, and only the limestone permeability is 
allowed to vary. If, however, the limestone is 
considered as of fixed permeability and the size of 
the fracture is varied, the effective resistance of 
the system follows a curve as I in Fig. 7. The 
limestone permeability has been taken here as 
0.01 darcy and the “reservoir radius’ as 300’. 

Curve II of Fig. 7 gives the fraction Q;/Q of 
the total flow, which enters the well directly 
through the fracture. 1—@Q;/Q, therefore, repre- 
sents the fraction which has been fed into the 
fracture by the surrounding limestone, within 
300’ of the well. It will be seen that for fracture 
widths greater than 0.75 mm, less than 14 percent 
of the total flow is fed into the fracture within 
300’ of the well. This explains the fact that for the 
large fractures the resistance of the composite 
system varies inversely as the cube of the fracture 
width, since in Eq. (1) the carrying capacity of 
the fracture itself, varies as the cube of its 
width, and when the fracture is large most of the 
flow passes into the fracture at distances from the 
well that are greater than 300’. 

Curve I of Fig. 7 contains the necessary infor- 
mation for computing the effect of acid treatment 
in limestone-fracture systems. For if we consider 
the effect of the acid to be essentially that of 
widening a fracture of initial width w;, by Aw, the 
ratio of the production capacity of the system 
after treatment, Q, to that before treatment, Qo, 
will simply equal the ratio of the ordinate of 
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Fic. 7. Resistivity of fracture-limestone systems for fixed limestone perme- 
ability (0.01 darcy). Curve I: pressure drop (atmos.) over first 300’ from well, for 
a unit flux (cc/sec.) per unit limestone thickness. Curve II: fraction of total flow, 
Q;/Q, entering the well which has traveled through the fracture for the whole 


distance of 300’. 


curve I at w; to that at w;+Aw.® The ratio 
Q/Qo, obtained in this manner, is plotted in 


* The assumption made in the analysis that the fractures 
are of infinite extent and are affected by the acid over their 
whole length should not cause large errors in the resultant 
values of Q/Qo, because of the high concentration about 
the well of the pressure distribution along the fractures. 
The exact value of the effective or real permeability at 


Fig. 8 as a function of Aw for various values of the 
initial fracture width w;. The limestone permea- 
bility has been taken as 0.01 darcy in all cases 
except that of the dotted curve, that for this case 
being 0.083 darcy. 


points distant from the well bore will, therefore, affect 
but little the over-all resistance of the system. 
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Fic. 8. Increase in production capacity of a fractured 
limestone due to acid treatment. Q/Qo= (production ca- 
pacity after treatment)/(production capacity before treat- 
ment). Aw =increase in fracture width. w;=initial fracture 
width. : limestone permeability =0.01 darcy. --—-: 
limestone permeability =1/12 darcy. 


Perhaps the most striking feature of these 
curves is the enormous effect possible, especially 
for the smaller fracture sizes, for even moderate 
increases in their widths.*® Although the values of 
Q/Qo of Fig. 8 must be somewhat high due to the 
assumptions’ of the infinite length of the frac- 
tures and their uniform widening by the acid, it 
seems certain, however, that even if strictly 
corrected for these assumptions, it would still be 
found that increases in production rate as high 
even as 100-fold can be reasonably explained by 
the fracture mechanism without invoking any 
ad hoc hypotheses. 

It will also be clear from Fig. 8 that the effects 
are larger for the systems whose fractures are 

* That increases in fracture width of 1 mm should not 
be uncommon follows from the fact that a single ‘‘dose”’ 
of 500 gallons of 18° Baume HC! can dissolve away a 
layer 4 mm thick from the side of an average limestone 
bed 10 ft. thick for a distance of 100 ft. 

1” The assumption of a single fracture will also tend to 
make the computed values of Q/Q» too high, for a given 


quantity of acid, but the orders of magnitude of the effects 
should remain the same. 
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initially narrower and hence were initially pro- 
ducing at lower rates. The fact that the dotted 
curve, which was computed for a limestone 
permeability of 0.083 darcy, is less than the 
corresponding curve (w;=0.1 mm) for a lime- 
stone permeability of 0.01 darcy, shows that if 
the initially low production rate is due to a low 
limestone permeability rather than small fracture 
widths the effect will again be relatively large as 
compared to that for a well with initially higher 
production rate (higher limestone permeability). 

In the case of the radial flow mechanism it was 
seen that the effects will be greatest for wells in 
which the inner zone permeability is a small frac- 
tion of that of the rest of the limestone, i.e., for 
those which show a relatively serious plugging. 
However, for a fixed inner zone permeability the 
wells with the highest permeability beyond this 
zone, and hence highest initial production rates, 
will show greater responses to the acid treatment, 
while if the permeability of the outer zone is 
considered as fixed, those in which the inner ring 
is of lowest permeability—and hence of lowest 
initial production rate—will show the greatest 
response. As field experience has shown that the 
smaller producers usually react best to acid 
treatment, it must be concluded that insofar as 
the mechanism may correspond to the radial flow 
case, these small producers are such because of 
serious plugging about the well bore rather than 
because of a low permeability in the main body of 
the limestone. 


SUMMARY 


The above results may be summarized as 
follows: 

1. Small increases in production rate—up to 
about 50 percent—due to acid treatment, may 
be explained on the assumption that the perme- 
ability of a small radial zone about the well bore 
has been increased from normal to higher values, 
as well as by removal of radial plugging or widen- 
ing of extended fractures fed laterally by the 
limestone proper. Unless the limestone does have 
extended fractures or is appreciably plugged near 
the well bore, acid treatment should be relatively 
ineffective in stimulating the production. 

2. Moderate increases in production rate—50 
to 500 percent—can be explained on a radial flow 
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basis only on the assumption that the wells were 
initially plugged, the extent of the plugging being 
the principal factor in determining the initial 
production rate, so that small producers will 
show larger responses. They can also be explained 
equally well by the assumption of the extended 
fracture flow. 

3. Increases in production rate appreciably 
larger than 500 percent, for wells of initially 
moderate production rates, can be explained only 
on the assumption that there are extended frac- 
tures in the limestone which are penetrated and 
widened by the acid. Here, the smaller producers 
should show the greater responses, whether their 
initially small production rates are due to low 
limestone permeabilities or small widths of the 
fractures. For very small producers, increases 
higher than 500 percent could also be explained 
on the radial flow mechanism of production, 
although it would have to be assumed that there 
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was initially a condition of almost complete 
plugging near the well bore. 

It should be pointed out that from preliminary 
“flow tests’’ one can make no a priori prediction 
as to the mechanism of production and hence 
foretell the probable effect of the acid. All 
mechanisms will show an approximately linear 
relation between the production rate and pressure 
differential from which only the resultant resist- 
ance can be derived; this resultant resistance can 
be synthesized from either the radial flow or 
fracture mechanism by an appropriate assign- 
ment of the many physical and geometrical 
constants available for defining the details of the 
system. 

The authors are indebted to Mr. M. Meres for 
assistance with the calculations of this paper, 
and to the Executives of the Gulf Research & 
Development Corporation for permission to 
publish it. 
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An exact theory is given of the seepage flux underneath 
dams of extended base width, with and without sheet 
piling, and under coffer dams with variable depths of 
excavation on the downstream side. The flux, per unit 
upstream-downstream potential drop, is expressed in all 
cases as 1/2 the ratio of two complete elliptic integrals of 
the first kind, the modulus of that in the numerator being 
the comodulus of that in the denominator, these moduli 
completely characterizing the detailed geometry of the 
various individual systems. In all cases the seepage 
becomes logarithmically infinite with increasing thickness 
of the underlying permeable stratum. The effectiveness of 
the sheet piling in cutting down the seepage flux is a maxi- 


mum when it is set at the heel or toe of the dam, and 
increases as the ratio of the thickness of the permeable 
stratum to the dam width increases. Although the flux de- 
creases with increasing piling depth, a large portion of the 
flux will in general persist until the piling penetrates to the 
very bottom of the permeable stratum. The same is true 
with regard to the penetration of coffer dams. The pressure 
distribution at the base of dams of finite width resting upon 
permeable strata is found to be only slightly different from 
that with infinitely thick underlying beds, except when the 
thickness of the permeable stratum is appreciably smaller 
than the width of the dam. 





INTRODUCTION 


EAVER' has given a rigorous treatment 

of the uplift pressures and moments acting 
on dams of extended base width, with and with- 
out sheet piling, resting on permeable strata of 
infinite thickness. Forchheimer? has calculated 
the seepage flux underneath dams of extended 
base width without sheet piling and underneath 
shut-off walls or coffer dams before any excava- 
tion has been made on the downstream side; in 
both these cases the finite thickness of the 
underlying permeable stratum was taken into 
account. 

Although, as will be seen below, the pressure 
distributions or uplift pressures found by Weaver 
will not be seriously affected by the finite 
thickness of the underlying permeable stratum, 
as long as the latter is not greatly exceeded by 


1 W. Weaver, J. Math. and Physics 11, 114 (1932). 

2 Ph. Forchheimer, Wien. Ber. 126, 409 (1917). An 
attempt to correct for the approximation in Forchheimer’s 
analysis in the case of the coffer dam problems (cf. footnote 
8) by adding to Forchheimer’s potential function a Fourier 
series potential, has been made by R. Hoffmann, Die 
Wasserwirtschaft 1, 174 (1934). The determination of the 
coefficients of this series, however, involves the laborious 
process of numerically solving sets of infinite linear equa- 
tions, which again can of course only be done approxi- 
mately. Hoffmann also applies his method to the problem 
of a dam with similar piling at both the heel and toe, the 
numerical work here being still more lengthy. While for 
practical purposes, Hoffmann’s results, and even Forch- 
heimer’s, are sufficiently accurate, the methods used here 
are not only exact, but they are also very much simpler 
than that of Hoffmann. 


the width of the base of the dam, and as long as 
the sheet pilings do not extend to the immediate 
vicinity of the base of the permeable stratum, 
the seepage flux implied by Weaver's analysis is 
infinite in every case. Forchheimer’s theory, on 
the other hand, which is based on a method of 
images, is not exact, even though it does give a 
remarkably close approximation to the correct 
values of the seepage flux. Furthermore, it is 
limited to the two special cases mentioned above. 
As a knowledge of the seepage flux underneath 
dams is of value in judging the dangers of 
undermining and underwashing of the dams and 
also, in the case of coffer dams, in estimating 
the capacity of the pumping equipment to be 
installed during the construction work, a more 
complete analysis of the whole problem was 
undertaken. In the following are given the results 
of an exact theory of the seepage flux, both 
under dams of extended base width and under 
coffer dams or simple sheet water shut-offs. 

As in the case of the previous work cited 
above, the dams will be taken as of great length, 
so as to make the problems two-dimensional. 
The analysis will be based upon the well-known 
method of conjugate functions, with the finite 
thickness of the permeable stratum leading to 
transformations involving the elliptic functions 
rather than the algebraic functions used by 
Weaver. The basis of these potential theory 
methods, on the other hand, is the fact that the 
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Fic. 1. Diagrammatic ‘‘z plane” representation of a dam 
without sheet piling, resting on a permeable bed of thick- 
ness H. 


flow of homogeneous fluids through porous 
media obeys, within such ranges of fluid velocity 
as occur in most practical problems, Darcy’s 
law, namely: 


t= —VO; b= (k/u)(p—vyegy), (1) 


where 0 is the macroscopic vector fluid velocity, 
@ is the velocity potential, k the permeability® 
of the medium, and is perhaps best defined by 
Eq. (1) itself, «4 and y. are the viscosity and 
density of the fluid, p the fluid pressure, g the 
acceleration of gravity, and y the vertical co- 
ordinate directed downwards. In the case the 
fluid is a liquid, an application of the equation 
of continuity gives at once, for homogeneous 
media : 

V’*b=0, (2) 


which reduces the problem to one in potential 
theory. 


I. DAMS OF EXTENDED BASE WIDTH. 
No SHEET PILING 


Representing diagrammatically a dam of ex- 
tended base width without sheet piling as in 
Fig. 1, an application of the Schwarz-Christoffel 
theorem gives for the function ¢(z), mapping the 
“s plane”’ figure EABCDE onto the upper half 
“t plane,’’ the relation 


s=c1f (dt/(t?—d*)) +c, (3) 


where the points E and D are to be mapped 


* For a complete discussion of the meaning of this term 
and the measurement of its numerical value, cf. R. D. 
Wyckoff, H. G. Botset, M. Muskat and D. W. Reed, 
Rev. Sci. Inst. 4, 394 (1933). 
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Fic. 2. ¢ plane representation of Fig. 1. 


% plane 














0 o, e" 





Fic. 3. ¢(#, WV) plane representation of Fig. 1. 


onto t= +d. Imposing the further correspondence 


z=0:-/=0; s=+4+w/2-t=+1. (4) 
where w is the width of the dam (base), and 
evaluating the constants in Eq. (1) and in- 
verting, it is readily found that 


tanh r2/2H 


j= , d=ctnh rw/4H, (5) 
tanh rw/4H 


H being the thickness of the underlying per- 
meable stratum. The resulting ¢ plane, with the 
boundary values of the potential and stream 
functions ®, V, is shown in Fig. 2. 

Now these boundary values of #, WV, give a 
¢=+7V plane diagram as in Fig. 3.4 Hence if 
this ¢ plane is mapped onto the ¢ plane of Fig. 2, 
we shall have the relation between ¢ and 2, 
through the intermediate variable ¢, giving the 
desired potential and streamline distribution in 
the original system. Applying once more the 
Schwarz-Christoffel theorem to effect this map- 
ping, we have 


dl 
eiteontanae +C2. (6) 


¢=84+71V=c, [-- — 
J ((#2—1)(#—d?))! 


Setting k=1/d, co=0, and choosing the limits of 


* The value of WV, is shown as negative as a consequence 
. the relation dV/dy = 0@/dx <0, for the coordinate system 
of Fig. 1. 





MORRIS MUS 











KAT 








Fic. 4. Seepage flux, per unit potential drop, (Q/A®), under a dam with no sheet piling as a 
function of the ratio of the dam width (w) to the thickness of the underlying permeable 


stratum (#7). 


the integral so that ¢=0-—B’, and noting that 
C”’ (¢=—%,)t=+1, c is found to have the 
value —%,d/K, where K is the complete elliptic 
integral of the first kind with modulus k.°® 
Hence, as the total potential drop in the system, 
A®, is 2%,, Eq. (6) can be finally rewritten as 


—A® ' dl 
f=o+iv= =a 
((1—f)(1—FF)) 
k=tanh rw/4H. (7) 








The seepage flux underneath the dam is 
evidently given by the value of %. This is 
given by the value of ¢ for D”’, corresponding to 
t=d=1/k. Thus: 


—Aeb pik dt 
—#,4i¥,=—— f 
2K (1—-f)A- 


— Ad 
=-——(K+1K’),* (8) 
2K 





—F#))! 


where K’ is the complete elliptic integral of the 


* Although the conventional symbol for the modulus of 
elliptic functions is k, the notation & will be used here in 
order to avoid confusion with the permeability constant 
(cf. Eq. (1)). 

* Whittaker and Watson, Modern Analysis, p. 502. 


first kind with modulus k’=(1—&?)!. As —¥, 
is evidently the total flux, Q, underneath the 


dam (per unit of dam length), it follows from 
Eq. (8) that 


RK'Ap kygK’Ah  _ 
= = ; k=tanh ww/4H, (9) 
2uK 2uK 





where Ap, Ah are the pressure and fluid head 
differentials corresponding to A®. 

The seepage flux per unit potential drop, i.e., 
Q/A&=K’'/2K,,’ is plotted in Fig. 4 as a function 
of w/H, the ratio of the dam width to the 
thickness of the permeable stratum. It will be’ 
seen that, as is to be expected, the flux, or 
effective conductivity of the system, decreases 
with increasing w/H from infinitely large values 
for w/H~0 to vanishing values for w/H-~. 
In an actual calculation of the seepage, one 
must, of course, multiply the ordinates of Fig. 4 
by kygAh/y.* 


7 Very complete tables of K, K’ and K’/K may be 
found in Tafeln der Besselschen, Theta, Kugel- und anderer 
Funktionen by K. Hayashi (1930). 

8 Forchheimer’s results which give values of Q/A® almost 
identical with those of Fig. 4, in spite of the approximate 
character of his analysis—his potentials are not strictly 
uniform along the y axis, as must be required by sym- 
metry—are: 
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_Fic. 5. Pressure distribution at the base of a dam with no sheet piling, for unit total pressure 
differential between upstream_and downstream sides. x/w=([distance from upstream side 
(heel) ]/[width of base of dam]. Curve I: Thickness of underlying permeable stratum (H) is 


infinite; Curve II: w/H=5; circles: w/H =1. 


Although Eq. (9) or its equivalent form, 
Q/Ab=K’'/2K, (10) 


has been derived here for the special case of the 
seepage underneath a dam with no sheet piling, 
it is to be observed that it is a perfectly general 
result for any system bounded by two segments 
of uniform potential alternating with two stream- 
line segments. The whole characterization of the 
particular system of interest with respect to 
the numerical value of Q/A® lies in the value of 
the modulus & of the elliptic integrals. For 
supposing the z plane representation of the 
original system has been mapped onto the upper 
half ¢ plane, this ¢ plane can be mapped onto a 
rectangle such as in Fig. 3 by an elliptic function 
transformation. If the latter is reduced to its 
normal form as in Eq. (7), the sides of the 
rectangle will be of length 2GK, corresponding 
to A®, and GK’ corresponding to AW, where G 
is a constant resulting from the normalization. 
Eq. (10) then follows at once. The modulus of 
' the elliptic integrals is given by the equation® 


k*(q— p)(r—s) +k | (g—p)?+(r—-s)? 


—(p+q—r—s)*}+(q—p)(r—s)=0, (11) 





1 cosh rw /8H—1 
i ae 
C/ae@ 2x log 3 cosh rw/8H—1’ 


for w/H <1, and Q/A@=1/[0.86+w/H], for w/H>1. 
a H. Bateman, Partial Differential Equations (1932), 
p. ; 





where p, q, 7, s, are the points on the real axis 
of the ¢ plane onto which are mapped the 
dividing points between the segments of uniform 
® and ¥ in the z plane. 

It should be further noted that Eq. (9) is not 
only applicable to the seepage underneath a 
dam without piling, but it should also afford a 
good approximation to the seepage around the 
banks of the dam if the latter are not directly 
anchored to impervious rocks. For the part of 
the dam above the downstream water level the 
pressure drop to be used in Eq. (9) may be 
approximated by the average of that on the 
upstream side to the depth of the downstream 
water level. The lateral seepage below the level 
of the downstream side will correspond to the 
total fluid head difference between the two faces 
of the dam. 

Although the main purpose of the above 
analysis has been that of deriving an expression 
for the seepage flux—Eq. (9)—it is also of 
interest to compare for this simple case the 
pressure distribution underneath the dam given 
here by the exact theory, where the finite 
thickness of the permeable underlying stratum 
is taken into account, and that given by Weaver 
for the case of an infinitely extended permeable 
stratum. For this purpose it is only necessary to 
notice that since ¥=0 along the base of the 
dam, Eq. (7) gives 





MORRIS MUSKAT 


z plane 

















y 


Fic. 6. Diagrammatic ‘‘z plane’’ representation of a dam 
with sheet piling, set to depth / in the permeable bed of 
thickness H. 


p ® 1 os 
—=— = ——-F(sin“ f," k); 
Ap Ab 2k 


tanh rx/2H _ 


a tanh rw, 4H 


t k=tanh ww/4H, (12) 


where ® and x are measured from the center of 
the dam. The resultant pressure distribution for 
w/H=5, with x measured from the heel of the 
dam and Ap=1, is plotted in Fig. 5 as the dotted 
curve, that for w/H=1 being indicated by 
the circles. The solid curve gives the plot of 
p=(1/7) cos (2x—w)/w, given by Weaver for 
an infinitely extended permeable stratum. It 
will be seen from these that the assumption of 
an infinitely thick permeable stratum underlying 
the dam gives a good approximation to the 
pressure distribution (and hence total uplift) 
along the base of the dam, even though it does 
imply an infinite seepage flux. 


II. SEEPAGE FLUx UNDER DAMs OF EXTENDED 
BasE WIpTH. SHEET PILING PRESENT 


When the dam is provided with sheet piling 
extending to the depth / in the permeable 
stratum of thickness H, the seepage flux can be 
again computed by the successive application of 
two conjugate function transformations derived 
from the Schwarz-Christoffel theorem. Thus, the 
z plane of Fig. 6 will be mapped onto the upper 
half of the ¢ plane with the correspondence 


A—t=-1; B-t=0; C-t=+1, (13) 


by the transformation function : 


¥:¥ E $-¢, F G $=, 0 ¥% 
-d -f “ ( 9g | 





Fic. 7. ¢ plane representation of Fig. 6. 


tdt 
s=c1f - +C2 
(—d?)(f—1)! 


Cy (d?—1)'—(f—1)! 
= log +¢2, (14) 
2(d?—1)} (d?—1)!+(f—1)} 








where the points D, E are supposed to go into 
t= -+d, and the sign of (#—1)! should be — for 
t<—1, + for ¢>+1, and + imaginary for 
—1<t<-+1. Imposing the correspondence of 
Eq. (13) and the requirement that at t=-d, 
z goes over into the line: z=x+7H, it is readily 
found that 


€¢:= —(2H/r) ctn rh/2H, co=0, 


d=cse rh/2H. (15) 


Eq. (14) can, therefore, be rewritten as 


?=1+ctn? (rh/2H) tanh? (rz/2H). (16) 


The ¢ plane diagram is indicated in Fig. 7 with 
the boundary values of ® and W, and the 
extremities of the dam, F’, G’ denoted by —f, g. 

To find the seepage flux, it is now only 
necessary to apply the results of the last section, 
which give 


Q/Ab=K’/2K, 

with d(g+f)k?—2(gf+d?)k+d(g+f)=0 

ef +d?—((d?—f?)(d?—g?))? 
d(g+f) ) 





hen 
where 


*=1+ctn? (rh/2H) tanh? (rw/2H)(1—Z/w); 
(18) 
f?=1+ctn? (rh/2H) tanh? (r%/2H), 


% being the distance of the piling from the heel 
of the dam. 

To show the effect of the position of the 
piling on the seepage flux, Q/A® is plotted in 
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Fic. 8. Seepage flux, per unit potential drop, (Q/A@), under a dam as a function of the 
position of the sheet piling. ¥/w= (distance of piling from heel)/(width of dam); (depth of 
piling) / (thickness of permeable stratum) = 1/2; (width of base of dam) / (thickness of permeable 


stratum) = 1. 



















































































Fic. 9. Seepage flux, per unit potential drop, (Q/A®), as a function of the depth of piling, 
(set at the center of the base). h//7=(depth of piling)/(thickness of permeable stratum); 
w/H =(width of dam)/(thickness of permeable stratum). 


Fig. 8 asa function of #/w for w/H=1, h/H=0.5. 
It will be seen that while the flux is a maximum 
for the symmetrical position of the piling, the 
total variation is not large. To see how Q/A® 
varies with the relative width of the dam, w/H, 
and piling depth, h/H, it will, therefore, suffice 
to suppose that the piling is set at the middle of 
the base. The values so obtained will furthermore 
give upper limits to the seepage flux, and hence 
be on the ‘‘safe side’ when used in calculations 
of practical flow systems. 

The variation of Q/A® with the depth of the 
piling, #/H, is shown in Fig. 9 for various values 


of dam width, w/H, and for #/w=1/2. As is to 
be expected, the seepage flux decreases with 
increasing piling depth, and dam width, although 
for dam widths appreciably larger than the 
thickness of the permeable stratum the flux is 
almost independent of the piling depth until 
the latter extends to the immediate neighborhood 
of the base of the permeable layer. The fact 
that the flux rises very sharply as h/H decreases 
from the value 1, shows that unless the piling is 
actually anchored to the impermeable base, a 
considerable fraction of the seepage without 
piling will still persist even when the piling 





MORRIS MUSKAT 


















































Fic. 10. Pressure drop over the piling (set at center of base), for a unit total drop across the 
dam, as a function of w/h, the (dam width) /(depth of piling). h/H (depth of piling) /(thickness 


of permeable stratum). 


penetrates 99 percent of the thickness of the 
permeable layer. 

While we shall not enter here into the details 
of the pressure distribution at the base of dams 
with sheet piling, as given by the analysis 
developed above, an indication of the accuracy 
of the more approximate theory of Weaver may 
be obtained by computing the pressure drop 
over the piling for the special case where the 
piling is set at the center of the base of the dam. 
By repeating the analysis leading to Eq. (7) it 
is readily found that this pressure drop, 6p, is 
given by 


1 al/g 
<a ((1—2)(1—R))! 
1 
=—F(sin='1/g,k), (19) 
K 


where 
k=e/d; g?=1+ctn* (rh/2H) tanh? (rw/4H); 
d=cse th/2H, (20) 


and Ap is the total pressure drop across the dam. 
Eq. (19) is plotted in Fig. 10 for various values 
of h/H. The curve for h/J7=0 corresponds to 
the case of an infinitely thick permeable stratum, 
where Eq. (19) reduces to 


bp 1 (w/2h)?*—1 


— =— cos~! ——____—_ 


Ap x (w/2h)?+1 


which may also be derived from Weaver's 
analysis. 

The fact that the deviations from this limiting 
curve do not become large until h/H>0.5, gives 
a justification for using the simpler theory for 
an infinitely extended permeable stratum in the 
study of the uplift pressures and moments in 
dams with sheet piling, at least for smaller 
values of 4/H. In any case this approximation 
should be ‘‘safe’’ from the point of view of the 
practical design of dams with sheet piling, as it 
will lead to uplift pressures and moments that 
are higher than those that will actually occur. 

As to the numerical values of the seepage flux 
that correspond to the curves of Figs. 4, 8, 9, 
as well as that presented in the following section, 
it should be observed that the numerical value 
of the rate of flow, Q, will be given by 


QO=(RLAP/p)(Q/A®), (22) 


where & is the permeability of the permeable 
stratum, uw the viscosity of the water, L the 
length of the dam, and Ap the pressure differ- 
ential between the upstream and downstream 
sides of the dam. If & is in darcys, L in cm, Ap in 
atmos., uw in centipoises, Q will be expressed in 
cc/sec. Thus, taking 


k=10 darcys; u=1 centipoise ; 
L=3048 cm=100’; 
Ap =33.9 ft. of HO =1 atmos., 
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Fic. 11. Diagrammatic ‘‘z plane” representation of a 
coffer dam set toa depth hy, in a sand of thickness H, with 
a downstream excavation of depth J. 


a value of Q/A®=0.1 corresponds to 
Q=3048 cc/sec.~48.32 gal./min. (23) 


In view of the magnitude of this seepage rate, 
and the fact that as shown by Fig. 9 values of 
Q/A®~0.1 persist for piling depths as great as 
99 percent of the thickness of the permeable 
stratum, it becomes clear that real safety from 
high seepage velocities and the dangers of 
undermining can only be attained, in the case 
of dams which are not set themselves directly 
on an impermeable base, by actually anchoring 
the piling to the impermeable bed. Unless the 
piling is so anchored it will cut down the seepage 
but little, except in such cases where the width 
of the dam is small compared to the thickness of 
the permeable stratum. 
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Ill. THE SEEPAGE FLUX UNDERNEATH 
CoFFER DAmMs 


While the results of the last sections should 
suffice to give the salient features of the seepage 
flux of water underneath permanent dam struc- 
tures, the use of temporary water shut-offs, as 
coffer dams, during the course of building 
hydraulic works lends practical interest also to 
the question of the seepage flux underneath 
such temporary devices for the shutting off or the 
reduction of the quantity of seepage water. A 
typical diagrammatic representation of such a 
coffer dam structure is indicated in Fig. 11, 
where the dam is set to the depth /, into the 
permeable stratum of thickness H/, the down- 
stream side having been excavated to the depth 
he, the’dam itself being taken as a thin imper- 
meable sheet—AB. Pumps are installed on the 
downstream side to remove the water seeping 
into the excavation. A knowledge of the seepage 
to be expected will be of value not only in 
estimating the proper capacity of the pumping 
installation but also in anticipating the dan- 
gers of underwashing and undermining the dam 
itself. 

Applying once more the Schwarz-Christoffel 
theorem, it is readily found that the z plane 
diagram of Fig. 11 will be mapped onto the 
upper half ¢ plane by the transformation function 
defined by 


((dy—1)(t+1))!+((d, +1)(¢—1))} 





dt C1 
wn 
(t—d;)(t+ds)(2—1)3 di+d, 


ds 
/. 
where the correspondence 


A-—-1; B-0; 


(d—1)' 
(d?—1)} 


C>+1; 





g 
((d,—1)(t+1))?—((d,+1) \(t—1))4 
((do+1)(t+1))!+((d ———) 
Ca, 
((do PE a ee 


D-d,; 


t——do, 


has been presupposed as indicated in Fig. 12. These conditions give successively : 


—d d» 
+ )+e. 
(d— 1)} 


TIC 
me (| 
d\+de 
(di—1)* +a 


Cy —d, 
=>— ~—| log 
(d,;+de) (d;7—1)} 





—_ 1)} 


the=Co. 


d,+1)} d> 
(d:—1)!—i(d+1)! 





ae) 
g — 
(d2+1)!—i(d2—1)} 


=o lo 
(d>—1)} 
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Fic. 12. ¢ plane representation of Fig. 11. 


Noting from these that c; is real, the condition 
that z=x+7H, for t>d,, gives, on equating the 
imaginary parts of both sides of Eq. (24) for 
t>d, : 


—7d\C 
H(d\+d:)=— 


+¢2(d,:+d2) /1. 
(d,7—1)} 


(27) 
By eliminating ¢c; and cz and introducing now 
the notation 


d;=seca;; d2=sec ao, 


(28) 


it is found that a; and a2 must satisfy the equa- 
tions 


sin ao/sin a;}=1—he/H, 


ai(1—he H)+ae=r7(1—hy, IT). 


(29) 


For the special case /2=0—no excavation on 
the downstream side of the dam—Eqs. (29) have 
the immediate solution: 


a;=ae=(r/2)(1—h,/H); 


d,;=d2=csc wh,/2H; ho=0. (30) 


In general, however, Eqs. (29) may be solved 
graphically without difficulty. 

The seepage flux may now be computed by 
means of Eq. (10), with & given by Eq. (11). 
By solving the latter, and using the notation of 
Eq. (28), it is found that 


cos (a@;+ae)/2 


Q/Ab=K'/2K; k= 


(31) 


cos (a@,;— az) /2 


For the special case where the downstream 
excavation has reached the depth of the piling, 


10 The solution for this case has also been obtained by 
Forchheimer (reference 8) to a high degree of approxima- 
tion—the error in the solution again involving a non- 
uniformity of the potential along By in Fig. 11—by a 
method similar to that he used for the problem of Sec. I. 
His final formulas are: 


Q/Ab=2/{4 log 2 ctn (#/4)(1—4,/H)}, for h,/H>1/2, 
and 


Q/Ab=(1/r) log 2ctn rh, /4H, for h,/H<1/2. 


h2=h,, the above solution breaks down, and the 
problem must be solved directly with the points 
C and B taken as coincident from the beginning, 
Carrying through such an analysis, in a manner 
exactly similar to the above, it is found that for 
this case k is given by 


k=1/(2H/h—-1); ho=hy. (32) 


The values of the seepage flux as given by 
Eqs. (30), (31) and (32) are shown plotted in 
Fig. 13 as a function of the fractional dam 
penetration, },/H, for various values of he/h,. 
Here, too, while the general decrease of Q/A® 
with increasing h,/H or decreasing he/h; is to be 
expected from physical considerations, one finds 
the marked persistence of the flux until the dam 
is set at the very bottom of the permeable layer. 
As the physical equivalent of the ordinates is 
here the same as that given by Eq. (22), the 
importance of actually anchoring the coffer dam 
to the bedrock will be evident. 

Although the values of Q/A® derived here and 
in the previous sections have referred to the 
total flux seeping into the sand along the whole 
infinite upstream bed EA and leaving along the 
whole infinite downstream surface CD, in Fig. 11, 
the finite limitations of these surfaces in the 
actual flow systems will not invalidate the 
analytical results. For it is not difficult to show 
that the very great part of the flux into or out 
of the porous medium is concentrated near the 
singular points A and C. Thus the seepage per 
unit length out of the surface CD will be given by 


dv dwv ;dx 


ana a 


- (33) 
dx dt 


Taking the case where /2=0, and d,;=d2.=d, 
it is clear that 


dy A dx Bt 


dt ((@-1)\(@—#))' dt (@—2)(2—-1)8 





Hence 
dv (d?—#*) 
—=const. ————- 
dx t 
ctn th/2H sech rx/2H 


(14+-ctn? rh/2H tanh? rx/2H)* 





(34) 


= const. 
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Fic. 13, Seepage flux, per unit total potentia! drop, (Q/A®), under a coffer dam as a function 


of the depth of the dam. h,/H=(depth of 


penetration of dam)/(thickness of permeable 


stratum). /t2/h, = (depth of downstream excavation) /(depth of penetration of dam). 


as the mapping of the z plane onto the ¢ plane is 


in this case effected by the same function as that 
defined by Eq. (16). It is thus seen that dW/dx 
rapidly approaches a vanishing value as one 
recedes from the dam by a distance of the order 
of or twice the thickness of the permeable 
stratum. In any case, however, the results given 
here are again on the “‘safe side’ from a practical 
point of view, in predicting values of the seepage 
flux that will exceed those in the physical system 
where EA and CD are actually finite. 


More complex problems, such as the seepage 
of water into caissons or under dams with two 
sheets of piling, can also be treated by the 
method used here provided elliptic integral 
transformations are used for the original mapping 
of the z plane diagrams onto the ¢ planes. A 
discussion of these problems will be given in a 
later paper. 

The writer is indebted to the executives of the 
Gulf Research & Development Corporation for 
the permission to publish this paper. 
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The Effect of an Electric Field on the Viscosity of Liquids 


E. D. ALcock, Guggenheim Aeronautics Laboratory, California Institute of Technology 
(Received November 8, 1935) 


Previous investigators had found that an electric field 
would change the viscosity of some liquids. This paper 
deals with this effect of an electric field on the viscosity and 
some of the factors upon which it is dependent. A rotating 
cylinder type of viscosimeter was used, the two cylinders 
being insulated from each other and the potential applied 
between them. Several liquids were investigated, some polar 


and some nonpolar. Nonpolar liquids showed zero effect, 
while polar liquids showed effects which were dependent on 
their polar moment. Changes in viscosity up to 200 percent 
were measured using fields up to 2000 volts/cm. It was also 
found that the effect was dependent on the velocity gradient 
present in the liquid. 





1. EXPERIMENTAL INVESTIGATION 


ECENTLY the effect of electric fields on the 
viscosity of liquids has been investigated by 
several men.' The results of these investigations 
seem to indicate that there is a change in the 
viscosity of some liquids when they are in an 
electric field. In general the method used by 
these men was to observe the quantity of liquid 
which would flow in a given time through a 
channel. The purpose of these early experiments 
was to determine whether or not there was an 
effect and which liquids showed this effect. No 
attempt was made to determine the dependence 
of the effect upon other variables such as field 
strength, stress in the liquid, etc. The results of 
these investigations seemed to show that in 
general the liquids with electric moments showed 
a slight increase in viscosity when subjected to 
the electric field that they used. Liquids with no 
electric moments were tried and showed no 
change. This paper deals with the investigation 
of the dependence of the effect upon field strength 
and stress in the liquid for several liquids. Since 
one of the factors to be investigated is the stress 
in the liquid, the viscosity should be measured in 
a flow where a uniform velocity gradient exists 
throughout the liquid. The type of viscosimeter 
chosen to carry out this work was a modified 
type of rotating concentric cylinders. The inner 
cylinder was insulated electrically from the outer 
cylinder and was charged to a high potential. The 
dimensions of the instrument can be seen in the 
scale drawing (Fig. 1). The velocity gradient in 
this type of instrument is not quite constant, the 
variation being about 10 percent with the dimen- 


' Herzog, Kudar and Paersch, Naturwiss. 36, 662 (1933); 
Herzog, Kudar and Paersch, Physik. Zeits. 35, 446 (1934). 


sions used in the present instrument. The guard 
cylinders were only twice as long as the gap 
across which the viscosity was being measured, 
the end effects, therefore, were considerable. The 
presence of these end effects disturbs the flow 
slightly and spoils the instrument for measuring 
absolute viscosities; relative viscosities, however, 
could be measured to within 2 percent. Therefore 
the data presented in this paper give the relative 
change in viscosity due to the electric field. The 
experiments were all carried out at Reynolds 
numbers which were below 2000, so that with the 
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Fic. 1. Viscosimeter. 
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Fic. 2. View showing suspension, lower guard cylinder and 
supports. 


exception of the secondary flow due to the end 
effects, the flow is laminar. The suspension con- 
sists of a section of the cylinder upon which is 
mounted a shaft carrying the mirror (Fig. 2). 
This is suspended by means of a steel wire from 
the top of the instrument and is free to rotate and 
is also free to swing away from the axis of the 
system. This type of suspension is unstable elec- 
trically so that the field strength which it was 
possible to use in this instrument was restricted 
in general to values below 2500 volts/cm. When 
the system becomes unstable and swings away 
from the axis, a bearing action takes place and 
the suspension tends to rotate in a direction op- 
posite to that of the outer cylinder. This tendency 
to rotate magnifies the effect on the readings due 
to the deviation from the axis and is at once de- 
tectable in the inability of the suspension to 
reach an equilibrium. Deviations of as little as 
0.001 inch from the axis cause an appreciable 
instability of the suspension. This fact necessi- 
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tates a very exact adjustment of the entire 
system so that the cylinders are coaxial. The 
suspension of the system has a moment of inertia 
of 2500 g cm? and the suspension had a torque 
constant of 1090 dynes/cm. When the liquid is in 
the instrument, it comes up inside of the center 
cylinder and furnishes a certain amount of 
damping to the suspension. The amount of 
damping furnished by the liquid is illustrated by 
the fact that without liquid the suspension is set 
into oscillation by stray air currents in the room 
and rarely ever is stationary while with liquid in 
the instrument, large oscillations are damped out 
in twenty to thirty vibrations. The gaps between 
the suspension and the guard cylinders are less 
than 0.5 mm, hence the torque transmitted 
through these gaps to the inner liquid is neg- 
ligible. 

The source of voltage in these experiments was 
a 6000-volt d.c. generator, the output of which 
was filtered before applying it to the instrument. 
This arrangement furnished a very steady source 
of voltage and could be easily regulated. Voltages 
were read by means of a voltmeter placed across 
the instrument and were accurate to about 20 
volts/cm. The deflections of the suspension were 
read by means of a telescope and a scale placed 
160 cm from the mirror. The scale readings are 
corrected to give the deflection in terms of cm of 
arc with a radius of 160 cm. The deflections of the 
suspension were restricted by the limitations of 
this method of observation to values less than 
30°. The accuracy of the readings is 0.5 cm since 
one cannot obtain an absolutely steady deflection 
due to minute variations in the speed of rotation 
of the outer cylinder. The outer cylinder is driven 
through a worm reduction gear by means of a 
shunt wound d.c. motor. The speed of the motor 
was regulated by means of external resistances in 
the field and shunt circuits. In this arrangement, 
the torque delivered by the motor decreases with 
the speed of the motor, hence the accuracy of the 
readings at low speeds is greatly decreased, vari- 
ations of more than 5 percent being observed 
over a period of two or three minutes. This 
limitation combined with the restriction on the 
deflections of the suspension, eliminated the pos- 
sibility of measuring the effect on liquids of vis- 
cosity higher than 0.04 poise with any great 
accuracy. 
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TABLE I. 
VISCOSITY DIPOLE MOMENT? 
LiguIp (Poises) (Debye units X 10"*) 
Carbon tetrachloride 0.0096 0 
Benzene .0065 0 
Chloroform .0056 1.15 
Ether .0023 1.14 
Nitrobenzene — 3.95 
Acetone .0033 2.73 
n-butvl bromide .0294 1.90 
Toluene .0058 40 
Xvlene-meta .0061 40 


No precautions were taken to control the tem- 
perature of the liquid, all measurements being 
made at room temperature. The time necessary 
for the measurement is of the order of one minute. 
Readings were taken before and after the applica- 
tion of the field, and in general showed no change 
so that the change in viscosity due to tempera- 
ture changes is negligible. The liquids used in 
these experiments (see Table I) were of ordinary 
commercial purity with the exception of the 
nitrobenzene which was specially prepared by 
Dr. A. N. Prater of the chemistry department of 
the Institute. 

The results of the measurements are given in 
Figs. 3a to 3g. The results of the experiments on 
the various liquids have been plotted graphically, 
the percentage increase in the viscosity being 
plotted against the voltage gradient across the 
liquid for various values of deflection of the 
suspension. The deflection of the suspension is a 
convenient measure of the stress in the liquid, 
since it is directly proportional to the shear 
transmitted by the liquid. The deflections men- 
tioned in this paper are centimeters of arc on the 
scale which was 160 cm from the mirror. Carbon 
tetrachloride and benzene showed no effect due 
to the application of the field. A reversal of the 
field gave no change in the effect. In general, the 
application of the field causes no change in the 
viscosity until a certain value of the field gradient 
is reached and then any further increase in the 
field causes an increase in the viscosity. In one 
or two cases, as in chloroform and n-butyl 
bromide at low values of the field the viscosity 
tends to decrease and then to increase. While the 
accuracy of the readings is not great in this 
region the author feels that there was a definite 


2? Taken from O. Fuchs and K. L. Wolf, ‘‘Dielectrische 
Polarization,” Hand und Jahrbuch der Chemischen Physik 
(1935). 
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_ Fic. 3. Experimentally determined percentage increase 
in viscosity with electric field. Ordinates are percent in- 
crease in viscosity. Abscissas are volts per cm. 


decrease in the viscosity in this region and it was 
not due to experimental error. Some of the curves 
show a tendency to reach a saturation in the 
increase of viscosity but the writer does not feel 
that the accuracy of the readings is sufficient to 
say whether this is real or just due to experi- 
mental error. In general the estimated error in 
any reading is about 10 percent and in many 
cases is much less. The increase in viscosity 
measured by the writer is much larger than the 
highest increase observed by other workers, and 
he feels that it is too large to be attributed to any 
other source than to the electric field. The trend 
of the results shows that the increase in viscosity 
for a given value of stress in the liquid is greater 
the larger the dipole moment of the molecule. 
The curves show that in every case as the viscous 
stresses increase, the percentage increase in the 
viscosity decreases for given strength of electric 
field and a given dipole moment. It is also inter- 
esting to note the similarity of the curves for 
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chloroform and diethyl ether, molecules which 
have the same dipole moments. The curves for a 
deflection of the suspension of 20 cm are identical 
for both liquids which would seem to indicate 
that the molecules had similar shapes as far as 
collisions with like molecules were concerned. In 
the case of toluene and xylene, however, the 
curves for the same deflection of the suspension 
are greatly different, even though the molecules 
have the same dipole moment. One might con- 
clude from this that the shape of the two mole- 
cules was different. 


2. A THEORY OF THE EFFECT OF AN ELECTRIC 
FIELD ON THE VIsCosITy OF LIQUIDS 


If the liquid is one which has a dipole moment, 
the electric field will tend to orient the molecules 
of the liquid in a direction which is parallel to the 
field. If E be the field strength and m the electric 
dipole moment of the molecule, the energy of the 
molecule due to the electric field will be Em cos 6, 
where @ is the angle between the dipole axis and 
the electric field. There are also thermal and 
viscous energies associated with the molecule. 
Thermal energy does not depend upon any angle, 
hence will cause no orientation. Viscous stresses, 
as shown by Raman and Krishnan in their theory 
of flow birefringence, tends to orient the mole- 
cules of the liquid at an angle of 45° to the direc- 
tion of the flow. This 45° orientation then is the 
position in which a minimum of shear is trans- 
mitted by the liquid, or in other words is the 
lowest value of the viscosity. The amount of ori- 
entation due to viscous stresses when in equi- 
librium with thermal energy is proportional to 
the product of the viscosity and the velocity 
gradient at the point, and is also dependent on 
the shape anisotropy of the molecule. If we now 
place the liquid flow in an electric field, we will 
then have an equilibrium between the three 
energies; thermal, viscous and electric. If the 
electric field is in a direction other than at 45° to 
the flow, the equilibrium position will be at some 
direction other than at 45° to the flow and will 
cause an increase in the shear which the liquid 
will transmit. Since the velocity gradient has not 
changed, this means that the viscosity coefficient 
has changed. 

Kuhn has also investigated the equilibrium of 
an elliptically shaped molecule in a viscous flow. 
He finds that the molecule is in rotation; the 
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angular velocity, however, is not constant but 
reaches a minimum when the major axis of the 
molecule makes an angle of 45° to the direction 
of flow. Statistically this is the same result that 
Raman and Krishnan obtain in a different 
manner. However, here we should notice that the 
molecule is in actual rotation, not merely oscil- 
lating about the 45° position as might be the case 
in Raman and Krishnan’s theory. If we now 
subject the molecule which has a dipole moment 
to an electric field, the rotation will cease when 
the field becomes sufficiently strong. The reason 
for this is that the dipole moment has a definite 
direction and the potential energy of the molecule 
has a minimum at only one angle while the poten- 
tial energy of the molecule with regard to the 
viscous forces has a minimum in two positions. 

The writer has calculated the torque due to the 
viscous flow upon a molecule making the assump- 
tion that the molecule is spherical and the result- 
ing torque is of the same order of magnitude as 
the torque due to the electric field at the value 
where the viscosity begins to increase. The cal- 
culation is only an approximate one and should 
be carried out by using a more exact formula for 
the torque due to viscous forces. 

Let us now see what predictions can be made 
from this simple qualitative theory as to the 
effect of an electric field upon the viscosity of a 
liquid. First we can say that if a molecule has no 
dipole moment the electric field should not affect 
the viscosity. This agrees with the results from 
the tests on carbon tetrachloride and benzene. 
Secondly we can say that the greater the dipole 
moment, the greater the effect, which is checked 
in general by the experimental results. Third, 
we can state that the greater the viscous force 
in the liquid, the less will be the effect of a given 
field on the viscosity, this is also checked by my 
experimental results in every case. Fourth, we 
can state that no appreciable increase in the vis- 
cosity can take place until the electric field has 
stopped the rotation. This also agrees very well 
with the experimental results, for in general the 
electric field causes little or no change in the 
viscosity until a certain critical field strength is 
reached. Above this critical point the viscosity 
increases rapidly as the field increases. 

In conclusion the author wishes to express his 
appreciation to Dr. von Karman of the Institute 
for his assistance and guidance in this research. 






















































2 —opeare mo 


MARCH, 





PHYSICS 


Design of Tuned Resistance-Capacity Coupled Amplifiers 








VOLUME 7 


E. A. Jounson, Department of Terrestrial Magnetism, Carnegie Institution of Washington 


(Received January 4, 1936) 


The behavior of resistance-capacity coupled amplifiers is studied with regard to their behavior 
as resonant networks. Expressions are obtained for the optimum value of circuit constants and 


the design procedure developed. Tuning at one cycle per second has been successfully applied to 


low noise amplifiers. 


T has been shown by D. G. C. Luck! that a 

resistance-capacity network of the type used 
in resistance-capacity coupled amplifiers is the 
equivalent of a parallel resonant LC-circuit. It 
is often important to design amplifiers for use 
at particular frequencies with a minimum band 
width to reduce the noise level and consequently 
the limit of measurement of the amplifier. The 
use of LC-circuits is often objectionable because 
of their bulkiness, feedback characteristics, in- 
stability, and sensitivity to external fields. The 
following analysis shows the limit of band width 
that can be obtained with RC-circuits that are 
not subject to these disadvantages. The circuit 
is shown by Fig. 1. ‘The equivalent circuit is 
shown by Fig. 2. 

Following Guillemin? 


/ —_ | 
= Aa) Zalar+ (Cart Ze) Zart Ze) (1) 


where Z,;=R,, Z2=Rz/(1+jRitCw), Z3=1/jCw, 
and Z,=R,/(1+jR,Cw). Then 


n/A =a+j(bw—d/w) 


Maximum amplification occurs at a frequency 
wo = (d/b)}, A max =pu/a = Ap, (4) 


0.707 of maximum amplification occurs at 
w= [(a?+4bd)!—a]/2b, we=((a?+4bd)'+a ]/26, 
and the band width is 


W =we—w,=a/b (5)e 


Or w= (wo?+(W/2)?)!— W/2, we= (woe? + (W/2)?)! 
+W/2), and the amplification is 


|A| =Ao/[1+(wo/W)%(x—-1/x)2}, (6) 


where x= w/w. 

These formulae are exact up to the frequencies 
where the real part of the input impedance of 
the tube becomes appreciable. Since this occurs 
at radiofrequencies, the following discussion is 
adequate for frequencies below 1000 kc. 

For the design of an amplifier with nearly 
constant amplification over a wide band, we 
observe that W/2 should be large compared to 
wo and that then 


we = W, w, = W(uo/W)?. (7) 



































|A| =p/(a?+(bw—d/w)? }}', (2) 
where To make w: large we make a large and 6 small; 
a=(1+R,/Ri+R,/R,+C,/C 
+R,C,/R,C+R,C,/R1C, Cn fF O— 
b=R,C(C,/C+C,/C+C,C,/C), =z on en 
e * pe, I= $y Ge; 
and d=[R,2/R.R,+R,/R,]/R,C. [ - 
Let r,=R,/Ri, 1,2=R,/Rz G=C,/C, and eo 
gg=C,/C. Then sie 
2, zZ 
A=1+ 8p +e t+ Got Wel a tol vs sd 
b= RyC(qv+dot+ re): (3) eg z z €g2 
d=r,(1+r,)/R,C. — oh : : “ 
~ FIG. 2 4 


1D. G. C. Luck, Proc. I. R. E. 20, 1401 (1932). 
?E. A. Guillemin, Electrical Engineering Department, 
Massachusetts Institute of Technology. 


130 


Fic. 1. Circuit of amplifier. 
Fic. 2. Equivalent circuit of amplifier. 
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RESISTANCE-CAPACITY COUPLED 


w. determines the upper frequency at which the 
amplification begins to drop off appreciably. 
The choice of tube is usually determined by the 
amplification desired per stage, noise level, etc., 
but it can be seen from (3) and (5) that a low 
product of R,(C,+C,) will be most favorable. 
It must be borne in mind that C, is the effective 
input capacity and is given approximately by 
C,= CoctMeCyg, when uw, =effective amplification. 


If R, is fixed, W is made large by making C 


large enough so that g,<r, and g,<r, and by 
making R, and R, small enough so that W is 
the desired value. w is adjusted by adjusting 
w) by means of C. In general R,(C,+C,) de- 
termines the merit of any particular tube for 
use in a broadly tuned RC-amplifier. The 
example given below is for a low noise amplifier, 
hence the choice of tube. 


EXAMPLE (1)—-TUBE OF TYPE 38 


E,=+9 volts; E,,=+6 volts; E.,=—1.5 
volts; J,=60X10-* ampere; J,=approximately 
10- ampere; C,=C,=10 mmfd; p=160; R, 
=1.210* ohm. 

Let R,=0.3X10® ohm, R,=1X10° ohm, 
C=0.4 mfd. Then r,=4; r,=1.2; qp=q,=2.5 
X10-°>;, a=14+44+1.2=6.2; b=1.2X10°X0.4 
x10-* 5K 10-°=2.410-*; and d=1.2(1+4)/ 
1.20.4=12.5. 

Also Ao =160/6.2=26; W=6.2/2.4x10° 
= 260,000 radians/sec. = 41,000 c.p.s.; wo = (12.5 
X10°/2.4)}=720 r.p.s.=115 c.p.s.; we=41,000 
c.p.s.; and w,;=0.324 c.p.s. 

The reduction in gain for four stages at 10,000 
c.p.s. is 1/[1+(10/41)? ?=0.89 and at 10 c.p.s. 
it is negligible. 

Instead of a wide band width, a narrow band 
. Is often desired to reduce the noise level or to 
reduce the effect of undesired frequencies present 
in the input. 

We will consider the three simplified cases 
where (I) g,=q,=q, (II) q,=0 and g,=g, (III) 
qp>=0 and g,=q. It is desired to find the condi- 
tions that will make W smallest. 


Case (I) 
From (3) and (5) 


* 


W=a/b=(1 +rp+r,)(1+q)/R,C(2q+¢q) (8) 
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and 


wo=(r,(1+7r,)/(2g+q") }*/R,C, (9) 
and we obtain 

W= (1 +r,+17,)/(r (1 +r») }} 

xl +q)/(2q+ 4")? Joo. 


Differentiating with respect to r, and setting to 
zero we get 


ro=1+7,. -_ 


— 


(10) 


By inspection g should be as large/as possible. 
If the value of [(1+q)/(2g+q@)*] be plotted 
against g it is seen that it is sufficient if 


q=0.3 
and we have W> 3w». 


(11) 


Z 

Case (II) = 
Similarly for Case (I1) g,=0 and g,=g. 
W=([(1+rp+1,)/Lro(1 +ry) Fai +roqht +r,)* LIs 


If we differentiate with respect to g and equate 
to zero .- : 


q=(1 +1rpt+1,)/T. 
W= (20 (1 +1 p+19)/(1+ 7p) Po. 
So the condition is that 2c. We 2:43 We 
ro<(1+17,), Wo rw», (13) 
wo =rg(1+r,)?/R,C(i+r,+r,)'. (14) 
Case (III) 
For Case (III) g,=0 and q,=q. We have 
from (3) and (5) 
W=[(1+r,+10)/(Lro(1 +r») }@') 
+(1+1p)'qt/rg* Joo. 


Differentiating with respect to g and equating 
to zero 


(12) 


Then 


q=(1+r,+r7,)/(A+r,) » > 
W=(2(1+1r,+1,)!/r,)wo- 


(15) 
(16) 


Thus we need to make (1+7,)<r, for a narrow 
band width and then W=2a . The condition 
that (1+7,)<r, is usually undesirable. This is 
because r, cannot be made very small because 
of excessive drop in R, if Rx, is made very large. 
If r, is made large by making R, small, amplifi- 
cation is sacrificed. Thus either Case (I) or 
Case (II) is preferable. 


and 
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TABLE I. 








RELATIVE 
Norse LEVEL 
(2 STAGES) 


(mfd) fo—hi| 





0.150 
015 
.0015 
.00060 
.00030 
.00015 
.000075 








In practice it is not usually possible to use 
Case (I) or Case (II) at will, although the latter 
is preferable. Since wo depends upon C which 
determines g, and g,, at very high frequencies 
it is not always possible to make ¢,<q,. In 
any case the minimum band width depends upon 
the resonant frequency and is given for Case (II) 
by fe—fi| =2fo since w=2rf. 

The design then is carried out as follows by 
using Case (II). Choose R, and R, observing 
(13), g is then determined from (12), and C is 
obtained from the desired value of fy by using 
(14). Finally C, is determined from the value of 
C and gq, since C,=Cq,. 

The values for C and C, of Table I are ob- 
tained by using Case (II) with the values 
R,=1.2X 10°, Rp =0.3 X 10°, Rs =1X 10°, g=6.2/ 
1.2=5.15, and Ay=10. 

It is evident that the desired sharpness can 
be obtained only at the expense of gain and is 
approximately one-half for a sharply tuned stage 
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3. 3. Frequency characteristic of five-stage amplifier 
tuned to one cycle per sec. 


compared to a broadly tuned stage with no 
shunting capacity. 

Fig. 3 shows the measured frequency char- 
acteristic of five stages of the amplifier calculated 
above, tuned to one cycle per second. This 
amplifier has a low noise level since noise level 
is proportional to square root of band width and 
this design has been applied to the measurement 
of the voltage of an alternating-current thermo- 
couple,* of Boltzmann’s constant from resistance 
noise,‘ and recently in a sensitive alternating- 
current photoelectric amplifier. 


3 Louis Harris, Phys. Rev. 45, 635 (1934). 
*C. Neitzert, Physics 5, 292 (1934). 





